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Résumé 


Dans cette thèse, nous avons mené une étude approfondie des propriétés physiques des 
matériaux bidimensionnels de la famille de trihalogénures de métaux de transition (MX3) dopées 
par des métaux de transition. L'approche de recherche utilisée repose sur la méthode des ondes 
planes augmentées (PAW), mise en œuvre dans le code VASP, dans le cadre de la théorie de la 
fonctionnelle de la densité (DET), en utilisant les approximations GGA et GGA + U pour 


prendre en compte les effets de corrélation. 


ES 


La première section de notre étude est consacrée à l'analyse des propriétés de la 
monocouche Mol; dopée avec différents métaux de transition, à savoir Sc, Ti, V, Cr et Mn. Les 
résultats obtenus montrent que le dopage avec le Sc et le Ti transforme la monocouche Mol; d'un 
semi-conducteur bipolaire ferromagnétique en un semi-conducteur ferromagnétique. En 
revanche, le dopage avec le V et le Cr donne des propriétés half semi-conductrices. Le dopage 
avec le Mn conduit à un caractère demi-métallique avec une magnétisation renforcée, même en 
présence du couplage spin-orbit (SOC). Enfin, tous les systèmes dopés présentent un 
ferromagnétisme robuste, en particulier la monocouche Mol; dopée au V, avec des températures 


de Curie élevées dans l'approximation du champ moyen. 


Dans la deuxième section de l'étude, les monocouches VCI; elles aussi sont dopées avec 
des métaux de transition et ont révélé des propriétés intéressantes. Le dopage avec le Sc et le Ti 
a produit des semi-conducteurs ferromagnétiques avec des gaps d'énergie indirects, tandis que le 
dopage avec le Cr a donné un semi-conducteur ferromagnétique avec un gap d'énergie directe et 
le dopage avec le Mn et le Fe a généré des semi-conducteurs sans gap de spin (SGS) et des semi- 
conducteurs ferromagnétiques bipolaires (BFMS), respectivement, suscitant un vif intérêt en 
spintronique. Tous les systèmes dopés ont affiché d'importants moments magnétiques, 
notamment 4,75 ug par cellule unitaire pour la monocouche dopée au Fe et une stabilité 
ferromagnétique accrue. En conséquence, ces monocouches de trihalogénures de métaux de 
transition dopées sont prometteuses pour les applications en spintronique en raison de leurs 


propriétés uniques. 


Mots clés: Matériaux 2D, Spintronique, Propriétés électroniques, Propriétés magnétiques, 


Calculs de premiers principes, Dopage. 


Abstract 


In this thesis, we conducted an in-depth study of the physical properties of two- 
dimensional materials from the family of transition metal trihalides (MX3) doped with transition 
metals. The research approach used is based on the augmented plane wave (PAW) method, 
implemented in the VASP code, within the framework of density functional theory (DET), using 


the GGA and GGA + U approximations to take into account correlation effects. 


The first section of our study 1s devoted to the analysis of the properties of the Mol 
monolayer doped with different transition metals, namely Sc, Ti, V, Cr and Mn. The results 
obtained show that doping with Sc and Ti transforms the Mol; monolayer from a ferromagnetic 
bipolar semiconductor into a ferromagnetic semiconductor. On the other hand, doping with V 
and Cr gives semiconductor properties. Doping with Mn leads to a half-metallic character with 
enhanced magnetization, even in the presence of spin-orbit coupling (SOC). Finally, all doped 
systems exhibit robust ferromagnetism, especially V-doped Mol; monolayer, with high Curie 


temperatures in the mean-field approximation. 


In the second section of the study, VCI3 monolayers were also doped with transition metals 
and revealed interesting properties.  Doping with Sc and Ti produced ferromagnetic 
semiconductors with indirect energy gaps, while doping with Cr produced a ferromagnetic 
semiconductor with a direct energy gap and doping with Mn and Fe generated spin-gapless 
semiconductors (SGS) and bipolar ferromagnetic semiconductors (BFMS), respectively, 
attracting Keen interest in spintronics. All doped systems displayed large magnetic moments, 
including 4,75 ug per unit cell for the Fe-doped monolayer and increased ferromagnetic 
stability. As a result, these doped transition metal trihalides monolayers are promising for 


spintronics applications due to their unique properties. 


Keywords: 2D Materials, Spintronics, Electronic Properties, Magnetic Properties, First Principles 


Calculations, Doping. 
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Introduction générale 


Durant les années 60, le domaine de la nanoscience n'existait pas encore malgré les 
progrès réalisés par la mécanique quantique concernant les principaux outils théoriques qui 
permettaient de comprendre le comportement de la matière à l'échelle microscopique. En 
1959, le célèbre physicien Richard Feynman, lors de sa conférence visionnaire intitulée « Il y 
a beaucoup de place en bas », mettait en lumière les possibilités passionnantes offertes par la 
réduction de la taille de la matière à l'échelle nanométrique et prédisait que les gens en l'an 
2000 se demanderaient pourquoi des efforts sérieux dans cette direction n'avaient pas 
commencé avant 1960 [1]. Bien qu'il soit difficile d'évaluer dans quelle mesure il a influencé 
la communauté scientifique, il est indéniable qu'en l'an 2000, il n'y avait plus aucune raison de 
s'étonner. Des progrès importants avaient été réalisés en suivant les idées de Feynman, jetant 
ainsi les bases de ce que nous connaissons aujourd'hui sous le nom de nanoscience et de 


nanotechnologie. 


Dans ses débuts, la nanoscience a jeté les bases du développement de diverses 
nanostructures, catégorisées en fonction du nombre de dimensions restantes au-dessus de 100 
nm : les nanoparticules (0 Dimension, OD), les nanofils (1 Dimension, 1D), les films minces 
(bidimensionnels, 2D) et les polycristaux massifs ainsi que les nanocomposites 
(tridimensionnels, 3D). Cette exploration des phénomènes quantiques à l'échelle 
nanométrique a conduit à des percées fondamentales dans les domaines de la nanoélectronique 
et de la spintronique, illustrées par la découverte de la Magnéto-Résistance Géante (GMR) 
[2,3] et de la Magnéto-Résistance par Effet Tunnel (TMR) [4] dans les films minces 
multicouches, révolutionnant ainsi les technologies de détection magnétique et de stockage de 


données. 


Ces avancées en nanoscience ont ouvert la voie au domaine florissant des matériaux en 
deux dimensions. Le graphène monocouche, en particulier, a déclenché une renaissance de la 
recherche grâce à ses propriétés électroniques extraordinaires et à son potentiel pour conduire 
la prochaine génération de dispositifs électroniques [5,6]. Au sein de cette sphère diversifiée, 
des matériaux en 2D, les aimants en 2D sont devenus récemment un point central de la 
recherche scientifique [7,8]. Il est intrigant de noter que le théorème de Mermin-Wagner 


prévoit théoriquement l'absence de magnétisation spontanée à longue portée dans les systèmes 
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bidimensionnels à des températures finies [9], cependant, ce théorème repose sur l'hypothèse 
du modèle de Heisenberg isotrope à deux dimensions. L'introduction de l'anisotropie 
magnétique au sein des matériaux bidimensionnels brise cette barrière théorique, ouvrant ainsi 


la voie à la réalisation de matériaux ferromagnétiques bidimensionnels. 


Dans le domaine des matériaux bidimensionnels, plusieurs familles de matériaux 
distinctes ont émergé. "Le graphène fonctionnalisé" englobe des matériaux tels que "le 
graphane", formé par une hydrogénation complète du graphène, et "le fluorographène", 
résultant d'une fluoruration totale, chacun présentant des propriétés uniques [10,11]. Une autre 
catégorie comprend les matériaux 2D du groupe IVA, tels que le silicène, le germanène et le 
stanène, composés respectivement de silicium (Si), de germanium (Ge) et d'étain (Sn), bien 
qu'ils soient confrontés à des défis en ce qui concerne leur stabilité sous forme autonome [12]. 
En revanche, la deuxième classe englobe les matériaux 2D dotés de structures en couches 
maintenues ensemble par de faibles forces de Van der Waals, notamment "le nitrure de bore 
hexagonal (h-BN)", connu pour sa stabilité, et "les dichalcogénures de métaux de transition 
(TMD)", offrant des propriétés électroniques réglables en monocouches et trouvant des 
applications en optoélectronique, en vallétronique et en spintronique [13,14]. Ces diverses 
familles de matériaux continuent de stimuler les avancées dans le domaine de la science des 
matériaux et de la physique de la matière condensée. D'autre part, la famille des 
trihalogénures de métaux de transition MX; (M : métaux de transition, X : atome halogène) 
[15,16] se démarque de manière significative grâce à son magnétisme intrinsèque [17-24]. En 
effet, l'une des premières réalisations expérimentales de matériaux ferromagnétiques 
intrinsèques en deux dimensions a été la monocouche de Crl:, qui est un membre de cette 
famille [25]. Cependant, la température de Curie de Crl:, qui est de 45 K, est bien inférieure à 
la température ambiante de 300 K, ce qui limite son application en spintronique [25]. Par 
conséquent, les chercheurs se sont orientés vers la modulation et l'amélioration des propriétés 
des membres de la famille MX; pour les applications spintroniques, car le matériau doit avoir 
certaines propriétés telles qu'une température de Curie élevée, un comportement semi- 


métallique, une polarisation élevée, etc [26-33]. 
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La réussite de la commercialisation des matériaux bidimensionnels dépend non 
seulement de leur efficacité, mais aussi de la capacité du laboratoire à produire de manière 
rentable des dispositifs électroniques et spintroniques à grande échelle. De plus, des études 
théoriques approfondies sont essentielles avant de considérer ces systèmes pour des 
applications pratiques, permettant ainsi une compréhension complète et une optimisation de 
leurs propriétés. Les calculs ab-initio représentent le fondement même de la physique 
computationnelle, car ils résolvent l'équation de Schrôdinger sans recourir à des données 
expérimentales, fournissant ainsi une prédiction d'une précision inégalée concernant la 
structure électronique et les propriétés des molécules et des matériaux. Ces calculs, basés sur 
la mécanique quantique, offrent des aperçus indispensables sur les géométries moléculaires, 
les énergies, les propriétés spectroscopiques et les mécanismes de réaction. Ils constituent des 
outils indispensables pour la compréhension approfondie et la conception de divers 
phénomènes chimiques et physiques, jouant ainsi un rôle central dans la recherche scientifique 
et le développement de nouveaux matériaux. Pour cela, nous avons utilisé des calculs de 
premiers principes sur la base de la théorie de la fonctionnelle de la densité pour étudier les 
propriétés structurales, électroniques et magnétiques des monocouches Mol: et VCI: dans le 


cas pur et dopé avec des métaux de transition 3d. 


Le manuscrit de cette thèse est organisé comme suit : 


Le premier chapitre propose une revue bibliographique approfondie des matériaux 
bidimensionnels dans un contexte plus large, mettant particulièrement l'accent sur la structure 
cristalline ainsi que sur les propriétés électroniques, magnétiques de la famille des 


trihalogénures de métaux de transition (MX3). 


Le deuxième chapitre aborde la base théorique de la théorie de la fonctionnelle de la 


densité, ainsi que les approximations employées dans cette étude. 


Le troisième et quatrième chapitre présentent les résultats des calculs concernant les 
propriétés structurales, électroniques et magnétiques des monocouches Mol, et VCI, à la fois 


dans leur forme pure et lorsqu'elles sont dopées avec des métaux de transition 3d. 
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Matériaux bidimensionnels : un aperçu 


L.1. Introduction 


Les matériaux bidimensionnels (2D) sont des matériaux cristallins constitués d'une seule 
couche d'atomes, où les électrons peuvent se déplacer librement dans le plan 2D, tandis que 
leur mouvement dans la troisième direction est limité et soumis à la mécanique quantique. 
Ces matériaux sont souvent classés comme des allotropes de différents éléments, 
principalement ceux qui possèdent deux liaisons covalentes [1]. Lorsqu'ils se présentent sous 
leur forme élémentaire ou composée, ces matériaux portent souvent dans leurs noms des 
suffixes tels que -ene, -ane ou -ide. Les matériaux 2D présentent des propriétés optiques 
fascinantes [2], une grande mobilité [3] et une large gamme de propriétés électroniques telles 
que la semi-métallicité [4], la semi-conductivité [5] et les propriétés supraconductrices [6]. 
Ces caractéristiques ont largement contribué à leur popularité croissante auprès de la 


communauté scientifique. 


L.2. Graphène 


Le graphène est l'un des matériaux 2D les plus étudiés en raison de sa grande popularité. 
Il se compose uniquement d'atomes de carbone et possède une structure simple. C’est un 
matériau allotrope du carbone, tout comme le diamant, le graphite et le fullerène et sa 


configuration électronique est [He] 2s°2p°. Il est constitué d'une seule couche atomique de 
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carbone arrangée en une structure en nid d'abeille. Cette structure a été découverte pour la 
première fois par Novoselov et Geim en utilisant une technique d'exfoliation mécanique à base 
de ruban adhésif, qui leur a valu le prix Nobel de physique en 2010. Cependant, cette méthode 
d'exfoliation mécanique est laborieuse et ne permet d'obtenir qu'une très faible quantité de 


graphène [7]. 


Depuis la découverte des premiers flocons de graphène, la production de cette substance 
s'est considérablement améliorée de manière exponentielle. Une alternative viable pour la 
synthèse de grandes surfaces de graphène est le dépôt chimique en phase vapeur (CVD) [8- 
10]. Cependant, la qualité du graphène obtenu par CVD était inférieure à celle obtenue par 
exfoliation mécanique [8-11]. Actuellement, la communauté scientifique s'efforce d'obtenir 


du graphène de haute qualité comparable à celui produit par exfoliation mécanique. 


Les électroniques basées sur le graphène offrent des perspectives prometteuses pour la 
recherche fondamentale et les applications pratiques et même des entreprises de renommée 
mondiale telles qu'Intel et IBM s'intéressent à cette technologie, car le silicium, qui a 
longtemps dominé le domaine de l'électronique, atteint ses limites et que le graphène s'avère 
être une alternative plus performante [12]. La mobilité élevée des porteurs de charge du 
graphène est l'un de ses principaux atouts pour une utilisation en électronique. Cette mobilité 
reste élevée même à des champs électriques élevés, ce qui est une caractéristique remarquable 
du graphène [12]. Cependant, il y a encore des défis à relever pour la communauté 
scientifique. Tout d'abord, la croissance de graphène de haute qualité à l'échelle industrielle 
reste difficile, malgré des progrès significatifs dans la croissance épitaxiale réussie [13,141]. 
Deuxièmement, les dispositifs électroniques basés sur le graphène doivent être capables de 
reproduire leurs propriétés pendant une longue période, afin de pouvoir servir le monde dans 
la vie quotidienne. Néanmoins, on peut s'attendre à de nombreuses applications du graphène 
dans un avenir proche, car l'utilisation du graphène dans les matériaux composites est une 
attente immédiate. Il convient de noter que la production de masse de poudre de graphène a 
déjà été réalisée dans une expérience, ce qui ouvre la voie à la préparation de matériaux 
composites à base de graphène à faible coût pour une variété d'applications [15]. Enfin, le 
graphène avec ses défauts a également suscité de l'intérêt dans une variété d'applications et 


avec de nombreuses propriétés intéressantes qui sont déjà connues et de nombreuses autres qui 
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sont attendues dans un avenir proche, il a le potentiel de révolutionner de nombreux domaines 


technologiques. 


Au fil des dernières décennies, l'intérêt s'est principalement porté sur le graphène, en 
raison de ses propriétés électroniques, optiques et mécaniques très prometteuses [12,16], 
cependant, l'absence de bande interdite dans le graphène limite son application dans de 


nombreux cas pratiques tels que la production de transistors [17]. 
L.2.1. Structure du graphène 


Le graphène possède une structure en nid d'abeille constituée d'atomes de carbone avec 
une symétrie de groupe ponctuel D&. Cette structure présente une distance interatomique plus 
courte (1,42 À) que celle du diamant (1,54 À), grâce à la présence de liaisons sp” dans le 
graphène au lieu de liaisons sp° dans le diamant. La cellule unitaire du graphène est constituée 
de deux atomes de carbone, notés A et B, tels que présentés dans la figure I.1, sachant que & 
représente la distance interatomique et les vecteurs de réseau direct et réciproque sont définis 


respectivement comme suit: 


di FF nn" 

| 2 BZ 

| 

" D 
se 


Première zone de 
Brillouin 


a) Réseau de Bravais b) Réseau reciproque 


Figure L.1 : Structure cristalline du graphène [18]. 
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L.2.2. Structure de bandes du graphène 


Dans les années 1920, les scientifiques pensaient que la synthèse des matériaux 2D 
comme le graphène était impossible. (Cependant, il a été démontré que le graphène se 
comporte comme un semi-conducteur, avec une bande interdite qui disparaît au point K de la 


zone de Brillouin, indiquant sa nature de semi-conducteur à gap nul (Figure [2). 


mm 


er 
Fermi energy 


<Nalence DT 
à ; Conduction ban 


Valence band 


Figure L2 : Structure de bandes du graphène près du niveau de Fermi. Un diagramme 
schématique en (A) 2D, (B) 3D et (C) cône de Dirac de K et K”, qui correspondent au niveau de 
Fermi de (B) [19]. 
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Contrairement aux semi-conducteurs conventionnels, la bande interdite du graphène devient 
linéairement nulle et forme une forme conique, connue sous le nom de cône de Dirac. En fait, 
les électrons dans le graphène suivent l'équation relativiste de Dirac avec une masse nulle, 
contrairement au mouvement non relativiste des électrons dans la matière condensée gouverné 


par l'équation de Schrüdinger avec une masse normalisée. 


Lorsque des électrons de faible énergie dans le graphène interagissent avec le potentiel 
périodique du réseau de graphène, de nouvelles quasi-particules émergentes peuvent être 
parfaitement décrites par des formulations de Dirac unidimensionnelles et bidimensionnelles 
avec une vitesse effective d'environ 106 m/s. Ces quasi-particules nouvellement apparues sont 
connues sous le nom de fermions de Dirac sans masse. Par conséquent, le graphène offre un 
moyen d'étudier les interactions fondamentales de la matière sans nécessiter d'accélérateurs de 
particules. En raison de ces phénomènes, le graphène pourrait être un candidat potentiel pour 
développer divers dispositifs électroniques. Cependant, le manque de bande interdite dans le 


graphène limite son utilisation en optoélectronique. 


La modification de la bande interdite du graphène est considérée comme une solution 
appropriée pour son utilisation dans les dispositifs optoélectroniques. Cependant, cette 
modification peut s'avérer difficile. Par conséquent, les chercheurs se sont tournés vers 
d'autres matériaux semi-conducteurs 2D pour résoudre le problème du gap nul du graphène. 
Dans les sections suivantes, l’une des familles de ces matériaux 2D couramment utilisés est 


brièvement décrite. 


L3. Trihalogénures de métaux de transition 


Les trihalogénures de métaux de transition MX; (M : un métal de transition, X : un 
halogène) ont suscité un intérêt considérable tant du point de vue théorique qu'expérimental. 
Ces matériaux présentent une grande anisotropie ainsi qu'un ferromagnétisme intrinsèque à 
l'échelle nanométrique, ce qui les distingue clairement des autres systèmes à faible 
dimensionnalité. De plus, leur couplage spin-orbite est particulièrement important pour 


ajuster les spectres de bandes électroniques [20]. 
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Bulk Monolayer Brillioun Zone 


Figure L.3 : Structure cristalline et zone de Brillouin du MX; massif et monocouche. (a) Cellule 
unitaire rhomboédrique de type Bil; du MX3 massif. (b) Vue de dessus de la monocouche MX3, 
avec la cellule unitaire hexagonale marquée par la ligne pointillée. (c) et (d) différentes vues 
latérales de monocouches MX3. (e) Zones de Brillouin du MX3 massif (ligne continue) et 
monocouche (ligne pointillée), correspondant aux cellules unitaires illustrées en (a) et (b) [19]. 


En 2017, le premier aimant 2D, Crl:, a été synthétisé avec succès au laboratoire [21]. 
Des calculs théoriques ont confirmé l'ordre ferromagnétique dans Crlz, ce qui a suscité l'intérêt 
pour de nombreux autres matériaux magnétiques 2D, tels que CrBr;, CrCls, VF3, MnP; et 
NiCl: [22,19]. La plupart des trihalogénures de métaux de transition adoptent une structure en 
couches stabilisée par des forces de van der Waals (vdW) faibles. Cette faible liaison inter- 
couches facilite l'exfoliation d'échantillons minces, telles que des monocouches et des 
bicouches, à partir de leur équivalent massif [23,24]. Sous forme de monocouche, toutes les 
structures de MX; partagent la même configuration atomique, où les atomes métalliques 
centraux (M) sont pris en sandwich entre deux plans d'atomes de X, créant un réseau 
hexagonal avec une coordination octaédrique partagée sur les bords. Cela signifie qu'il y a six 
atomes de ligands (X) les plus proches voisins entourant les atomes de M centraux [19,23-25]. 
Afin de voir la structure optimale, la cellule unitaire des structures de monocouches MX; (X = 
F, CI, Br, D a été présentée dans la figure L.3, comprenant des vues de dessus et de côté. La 
cellule unitaire en phase massif et en monocouche est formée de deux unités, chacune 


contenant un atome M et trois atomes d'halogénure. 
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L.3.1. Synthèse des matériaux MX; en monocouches 


Dans le domaine des matériaux 2D, un aspect crucial réside dans la synthèse et la 
caractérisation de ces structures remarquables. Plusieurs méthodes expérimentales ont été 
développées pour la synthèse de ces matériaux, tels que le dépôt chimique en phase vapeur 
(CVD), l'épitaxie par jet moléculaire (MBE) et les méthodes basées sur des solutions [26-28]. 
Ces approches permettent une croissance contrôlée des monocouches tout en offrant une 
grande flexibilité dans le choix des précurseurs et des conditions de croissance, ouvrant ainsi 
la voie à leur exploration et à leur utilisation dans une variété d'applications [29,30]. Une fois 
synthétisés, ces matériaux nécessitent une caractérisation approfondie afin de comprendre 
leurs propriétés et leur comportement. Des techniques de caractérisation telles que la 
microscopie, la spectroscopie et l'analyse par diffraction fournissent des informations 
précieuses sur la structure, la morphologie et les propriétés électroniques de ces matériaux 
[31,32]. Cependant, avant de se plonger dans ces techniques, il est important de considérer la 
faisabilité expérimentale de la synthèse et de la manipulation de monocouches. La faisabilité 
de la synthèse de monocouches, à partir de leurs structures massives, est déterminée en 
évaluant leur énergie de clivage [23]. Cette énergie est définie comme l'énergie nécessaire 
pour séparer ou fendre un matériau le long de plans ou d'interfaces spécifiques. En analysant 
la différence d'énergie entre les états massif et fracturé à des séparations intercouches 
croissantes, on peut évaluer l'énergie de clivage. Le tableau I.1 présente quelques valeurs 
d'énergie de clivage de différentes monocouches MX3. Les valeurs d'énergie de clivage 


inférieures à celle du graphite (0,37 J/m°) suggèrent leur faisabilité expérimentale. 


En 2017, la synthèse des trihalogénures de métaux de transition a abouti à l'obtention de 
flocons de Crl; atomiquement minces par exfoliation mécanique sur substrats de silicium 
oxydé [21]. Après ça, la méthode de transport chimique en phase vapeur a permis la 
croissance de cristaux monocouches de VI:, mettant en évidence son potentiel pour de futures 
explorations et applications [34]. En s’appuyant sur ces avancées, en 2022, une percée a été 
réalisée avec la croissance directe de grandes écailles de CrCl; jusqu'à une épaisseur 
monomoléculaire grâce à la technique de transport physique en phase vapeur (PVT) [351]. 
Cette méthode a permis un contrôle précis de l'épaisseur des écailles, ouvrant ainsi des 


possibilités pour l'étude et l'exploitation des propriétés uniques des matériaux 
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antiferromagnétiques. Très récemment, Mastrippolito ef al. ont réussi l'exfoliation mécanique 
in-situ de flocons minces de VCI:, permettant ainsi l'étude de leurs propriétés magnétiques en 
couches [36]. Ces avancées significatives dans la synthèse et la caractérisation des 
monocouches des trihalogénures de métaux de transition ont ouvert des perspectives 
prometteuses pour la recherche future et les applications pratiques ce qui offre de nouvelles 


opportunités pour explorer les propriétés uniques de cette famille de matériaux. 


Tableau LI : Energie de clivage de quelques monocouches MX; [33]. 


Energie de clivage (j/m°) 


Matériau 
Rhomboédrique  Monoclinique 

CrCB 0,10 0,13 
CrBr; 0,14 0,14 

Cri 0,16 0,16 
VC OA - 
VBr; 0,28 - 

VB 0,18 - 
NiBrs 0,14 - 
FeBr; 0,12 - 


L.3.2. Propriétés électroniques des monocouches MX; 


La compréhension des propriétés électroniques est cruciale dans la découverte de 
nouveaux matériaux. Les rôles distincts des semi-conducteurs et des métaux dans les 
applications sont étroitement liés à leur structure électronique des bandes d’énergie. Cette 
structure, qui est l'un des aspects clés en science des matériaux, permet de décrire les 
propriétés électriques, optiques et magnétiques complexes. Les orbitales d partiellement 
remplies des ions de métaux de transition dans les monocouches MX; permettent des états 
d'oxydation variables ce qui conduit à des structures cristallines diverses. L'hybridation des 
orbitales d des métaux de transition avec les orbitales des halogènes conduit à la formation de 


haisons covalentes, ce qui leurs confère des propriétés métalliques ou semi-conductrices. 
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Des études récentes ont révélé que les trihalogénures de Cr, Fe et Ti présentent un 
comportement isolant, tandis que les trihalogénures de Ni, Mn et V affichent un état 
électronique unique appelé demi-métal de Dirac [33]. Ce dernier se caractérise par un 
comportement métallique avec la présence d'un cône de Dirac dans l'un des canaux de spin 
[38,39]. Un exemple notable est illustré dans la figure [.4 (a)-(b) pour les structures de bandes 
électroniques des monocouches VX3 (X = IL Cl). L'analyse de cette figure révèle que le canal 
de spin-up présente une importante bande interdite, tandis que le canal de spin-down affiche 
un cône de Dirac dans le diagramme de bandes. Ces structures présentent des caractéristiques 
demi-métallique de Dirac et ont un potentiel en tant que dispositifs de filtrage de spin [40,41]. 
De plus, le diagramme de bandes montre une relation de dispersion linéaire près du niveau de 
Fermi, ce qui confirme la présence de fermions de Dirac sans masse. Ces derniers jouent un 
rôle crucial dans l'émergence de phénomènes topologiques dans les matériaux. Ces 
phénomènes, tels que les isolants topologiques, les semi-métaux de Weyl et les semi-métaux 
de Dirac, présentent des propriétés uniques et prometteuses pour des applications en 


informatique quantique et en spintronique [42-441]. 


Euergy(eV) 


LEE 
RS Y 


1 M K r DOS(a.u.) r M K r DOS(a.u.) 


Figure L.4 : Structure de bandes électroniques polarisées en spin avec les densités d’états (DOS) 
des monocouches (a) VC; et (b) VI:[37]. 
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L.3.3. Propriétés magnétiques des monocouches MX; 


Les chercheurs ont toujours été fascinés par le magnétisme à basse dimension. En 
réduisant la dimension de trois à deux, les fluctuations thermiques augmentent 
considérablement. Avec la diminution de la dimension spatiale, il a été constaté que les 
phénomènes magnétiques deviennent dépendants de la dimensionnalité du spin [45,46]. Selon 
le théorème de Mermin-Wagner, à une température finie, l'aamantation spontanée est 
impossible à atteindre dans les systèmes 1isotropes à courte portée à deux dimensions [47]. 
Malgré l'induction réussie de magnétisme dans certains matériaux 2D non magnétiques grâce 
à des techniques telles que le dopage, l'ingénierie de contrainte, l'ingénierie des défauts et 
l'effet de proximité, le contrôle précis de ces méthodes en laboratoire reste un défi persistant. 


En conséquence, les chercheurs se sont tournés vers la recherche de nouveaux matériaux 2D 


qui présentent un magnétisme intrinsèque. 


Les prédictions théoriques de nombreux matériaux ferromagnétiques intrinsèques, tels 
que CrGeTe: [48], 1T-TaN: [20] et CrX3 (X = F, CI Br, D [21], ont suscité un intérêt 
considérable dans ce domaine. Parmi eux, la synthèse de la monocouche de Crl; [22] se 
distingue comme une réalisation notable. Cette découverte cruciale démontre non seulement 
l'existence d'un ferromagnétisme intrinsèque dans les systèmes 2D, mais ouvre également des 
voies d'exploration pour étudier plus en profondeur les propriétés magnétiques des 


trihalogénures de métaux de transition. 


La température de Curie Tc est l'une des caractéristiques les plus importantes des 
matériaux magnétiques. Elle est simplement définie comme étant la température au-dessus de 
laquelle certains matériaux perdent leurs propriétés magnétiques permanentes. Les valeurs de 
Tc, présentées dans le tableau [.2, sont estimées en utilisant l'approximation de la théorie du 
champ moyen, une approche qui permet d'estimer la température critique (Tc) des 
monocouches de MX3. Bien qu'elle puisse surestimer les valeurs réelles, cette méthode offre 


un moyen de comparaison entre les différentes monocouches de MX:3. 


De manière intéressante, nous observons que les MX; à base de Mn et de Ni possèdent 
des valeurs de Tc plus élevées par rapport aux autres monocouches MX3. Le tableau [2 
affiche également les valeurs de l'énergie d'anisotropie magnétique (MAE), qui est définie 
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comme étant l'énergie requise pour modifier l'orientation du moment magnétique dans un 
matériau, passant de son axe préférentiel, appelé axe facile, à une direction différente, appelée 


axe dur. La MAE a été calculée théoriquement et expérimentalement pour quelques structures 


Tableau L.2: Energie d’anisotropie magnétique (MAE) par cellule unitaire (*: MAE en 
erg/cm”, **: MAE en meV), température de Curie Te et moment par atome métallique de 
quelques monocouches MX; [33]. 


Matériau MAE par atome Te (K) Etat fondamental Moment 
métallique (HeV) magnétique magnétique (LB) 
CrF; 119,0 41 FM 29 
CrCB 24,68 29,7 FM 2,12 
CrBr; 159,54 41,3 FM 3,25 
Cr 803 46,4 FM 3,44 
MnF; -0,013* 450 FM 3,92 
MnCI -0,46* - FM 4,08 
MnBr; -8,71* - FM 4,18 
Monk -11,86* 720 FM 4,27 
NiCB 217 497 FM 1,15 
NiBr; 263 595 FM 1,18 
Nils -320,6 682 FM 1,23 
FeF; 0,005* 56 FM 4,35 
FeCl 0,055* 716 FM 4,10 
FeBr; 0252 116 FM 4,01 
Fel: 1,145* 148 FM 3,88 
VC - 80 FM 2 
VBr; 300 190 FM 1,97 
VB 0,29* 98 FM 2 
TiClz - 5 ; : 
TiBr; 0,8 15 FM 0,98 
Ti 297 90 FM 1,04 


monocouches MX; (Tableau I.2) et ces valeurs sont suffisamment élevées par rapport aux 
matériaux ferromagnétiques massifs tels que Fe (1,4 eV par atome) et Ni (2,7 eV par atome) 
[49], garantissant une meilleure efficacité dans le stockage des données magnétiques contre 
l'instabilité thermique et elles sont également comparables aux valeurs de 0,56 erg/cm” de 
l'hétérostructure Au/FeCo/M£eO pour les applications de mémoire vive magnétoélectrique 


(MeRAM) [50]. 
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L.3.4. Modification des propriétés électroniques et magnétiques des MX; 


Diverses méthodes et approches ont été utilisées pour ajuster et améliorer les propriétés 
électroniques et magnétiques des monocouches MX; en vue d'applications pratiques. Les 
structures de bandes des monocouches de Crl: ont été étudiées sous l'effet de contraintes 
uniaxiales et biaxiales de compression et de traction. Il a été démontré que Crl; dans son état 
initial est un isolant dans les deux canaux de spin-up et de spin-down. Sous une contrainte 
uniaxiale, la bande interdite du canal de spin-up se dégrade à 2 eV à 2 % de contrainte et reste 
pratiquement constante jusqu'à 12 % de contrainte, tant qu’en compression que qu'en traction. 
Dans le canal de spin-down, la bande interdite diminue sous une contrainte uniaxiale de 
compression et augmente légèrement sous une contrainte uniaxiale de traction. Ces variations 
des valeurs de la bande interdite s'expliquent par le chevauchement des orbitales lorsque les 
atomes se rapprochent ou s'éloignent. Les configurations asymétriques de contraintes 
biaxiales ont conduit à des bandes interdites réduites, tandis qu'un comportement oscillant a 
été observé dans une configuration pour le canal de spin-up. Dans l'ensemble, la manipulation 
de la contrainte a permis de contrôler significativement les propriétés électroniques de la 


monocouche de Crl: [51]. 


D'autre part, les expérimentations ont clairement démontré que la technique d'alliage est 
très efficace pour ajuster considérablement la bande interdite du composé monocouche 
Cle ,Br,. En effet, il a été remarqué que la bande interdite augmente proportionnellement 
avec l’augmentation de la concentration de Br, permettant ainsi sa variation de 0,86 eV à 1,70 
eV. Pour des concentrations de Br plus faibles (x < 2), les composés alliés se comportent 
comme des semi-conducteurs à bande interdite directe. En revanche, pour des concentrations 
plus élevées, une transition vers un semi-conducteur à bande interdite indirecte est observée 


(521: 


D'un point de vue structural, on peut en conclure que l'application d’une contrainte 
modifie la longueur et les angles de liaisons, ce qui affecte de manière positive la température 
critique de transition correspondante (T.). Il a été également observé qu'une augmentation de 
la contrainte améliore Tc, et qu’en appliquant une contrainte de traction biaxiale de 10%, il est 


possible d'atteindre une valeur de T< de 130K pour Crl; [24]. De plus, il était remarqué que 
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l'énergie d’anisotropie magnétique (MAE) de Crl; augmente sous contrainte de compression 


biaxiale, tandis qu'une tendance inverse a été observée pour CrCl et CrBr; [53]. 


De plus, le dopage avec les porteurs de charges est une autre méthode efficace pour 
moduler ou ajuster les propriétés magnétiques. On a observé que les paramètres d'interaction 
d'échange se sont améliorés par le dopage en trous (jusqu'à 0,75 trou par électron par maille 
unitaire) et le dopage en électrons (jusqu'à 0,75 électron par trou par maille unitaire) pour VCI; 
et VI, respectivement, ce qui a permis de démontrer une amélioration supplémentaire de la 
température critique de transition jusqu'à la température ambiante (353K et 346K pour VC et 
VB, respectivement) [37]. En revanche, l'effet combiné de l'ingénierie de contrainte et de 
l'alliage sur monocouches alliées Crl_,Br, a été étudié et le moment magnétique total trouvé 
n’était pas très sensible à la variation de la concentration x en Br et à l'application de 
contraintes de traction et de compression. De plus, Tc augmente légèrement avec la contrainte 


de traction [52]. 
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Chapitre IT : 
Méthode de calcul 


IL.1. Introduction 


L'électron et le noyau sont les constituants de base d’un matériau et sont responsables 
de sa nature physique et chimique ainsi que de ses différentes caractéristiques atomiques et 
moléculaires, telles que : les propriétés optiques, magnétiques, de transport et les structures 
cristallines qui dépendent essentiellement de la structure électronique. Par conséquent, l'étude 
des structures électroniques a toujours attiré l'attention de la communauté des sciences des 
matériaux, notamment en physique et en chimie. Cependant, les solutions de la structure 
électronique peuvent ne pas être si simples, car les interactions électroniques sont de nature de 
mécanique quantique dans la matière et leur description devient plus complexe en raison du 
nombre croissant d'électrons. Ce problème conduit à la physique à plusieurs corps, qui est 


considérée comme une branche importante de la physique des matériaux. 
IL.2. Problème à plusieurs corps 


Le système à plusieurs particules est exprimé par des fonctions d'onde tout-électron 
T,R, ,t), qui dépendent généralement de la position et du temps. L'équation (IL1) de 
v Ra q P 8 P P q 


Schrôüdinger contrôle la dynamique des systèmes non relativistes; 
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ih— = AY (IL.1) 


Ici, l’hamiltonien (H) représente l'opérateur d'énergie totale ayant la forme suivante 
pour le système à plusieurs corps en interaction, qui comprend un certain nombre d'électrons 


et de noyaux en interaction, 


PV D Dr DDR R| DD ce (AL2) 


où la masse de l'électron et du noyau est représentée par m, et M,, respectivement, et la 


me ième 


position respective du i% électron et du @ noyau est représentée par T7; et ha 
respectivement. Le numéro atomique du noyau correspondant est représenté par Z,. Dans 
l'équation (IL2), la première et la deuxième partie montrent l'énergie cinétique des électrons et 
des noyaux, respectivement. Le reste des termes indique l'énergie potentielle, issue des 
interactions entre électron-électron, électron-noyau, et noyau-noyau, respectivement. Il n'y a 
pas de terme explicite dépendant du temps dans l'hamiltonien ci-dessus. Par conséquent, la 
fonction d'onde peut être simplement écrite comme un produit des termes spatiaux et 
temporels, #({r, Ré}, t) = Er, Rien, ce qui conduit l'équation (IL.1) à une forme 


simple indépendante du temps : 
Ây = Ey (IL3) 


où E représente l'énergie totale du système. Cependant, la solution de l'équation de 
Schrôüdinger se limite à un nombre restreint de systèmes. Il est donc nécessaire d'intégrer des 


approximations pour appliquer l'équation de Schrôüdinger à tous les types de systèmes. 


L’approximation de Born-Oppenheimer [1] indique que le mouvement des noyaux est 


beaucoup plus lent que celui des électrons en raison de la nature des noyaux qui sont —10° 
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plus lourds que les électrons. Aüïnsi, il est probable que les électrons puissent adopter la 
position instantanée des noyaux très rapidement sur le temps de déplacement des noyaux. Par 
conséquent, les fonctions d'onde des noyaux seront indépendantes des coordonnées 
électroniques et se divisent en deux termes, à savoir le terme des noyaux et les termes 
électroniques. Ainsi, la partie séparée des noyaux et la partie électronique de l'hamiltonien est 


obtenue comme suit : 


2 


BDD) ED) 
= | — a | Tan II.4 
2Me : É _— | — Ra | _ R;| ( ) 


À j>i 


La méthode d'Ewald [2] est utilisée pour traiter l'interaction entre les noyaux, le dernier 
terme de l'équation (11.2) est ajouté à cette interaction pour estimer l'énergie totale du système. 
Même après avoir considéré l'approximation ci-dessus, la résolution de l'équation de 


Schrôdinger n'est pas une tâche facile pour les raisons suivantes : 


e Il ya -107 nombre d'électrons, donc un total de 4N variables est nécessaire pour 
estimer la fonction d'onde du système à plusieurs corps, 

e La présence d'électrons dans les solides, qui affectent le mouvement des autres 
électrons via le terme de corrélation (électron-électron), le troisième terme de 


l'équation (11.4). 


Par conséquent, pour trouver une solution possible, différentes approches ont été 
proposées pour résoudre le problème à plusieurs corps. Parmi lesquelles, on trouve celui de 
Hartree [3] qui a introduit la toute première approche pour résoudre le problème ci-dessus en 
considérant la fonction d'onde à plusieurs électrons comme un produit de fonctions d'onde à 
un seul électron. Les équations hamiltoniennes à une seule particule peuvent être estimées en 
utilisant le principe variationnel. Les équations obtenues sont les mêmes que celles de 
l'équation de Schrüdinger avec un terme supplémentaire de potentiel "Hartree”". Cependant, il 


manque la description de la fonction d'onde fermionique. 
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Hartree et Fock [3] ont résolu ce problème en considérant la fonction d'onde 
électronique à plusieurs corps sous la forme du déterminant de Slater. Des équations 
hamiltoniennes similaires sont obtenues en utilisant la méthode variationnelle. Dans ce 
formalisme, un nouveau terme, appelé potentiel d'échange, a été introduit avec le potentiel de 
Hartree. Cette approche est utilisée avec succès pour un système fini. Cependant, aucun effet 
de corrélation électronique n'a été pris en compte dans ce formalisme et il reste donc imprécis. 
Les formalismes mentionnés ci-dessus (Hartree et Fock) sont basés sur la fonction d'onde des 
électrons à plusieurs corps. Par conséquent, on s'attend à ce qu'ils soient très coûteux en 
termes de calcul pour les systèmes de grande taille. La fonction d'onde ne peut être mesurée 
expérimentalement en raison de sa nature complexe, car elle dépend de 4N variables ; trois 
variables spatiales et une variable de spin pour chaque N électrons. Le degré de liberté a été 
réduit par la densité électronique (s1 elle est utilisée comme variable), ce qui permet de réduire 


considérablement le temps de calcul. 


Résoudre l'équation de Schrüdinger pour un système à plusieurs corps en utilisant la 
densité électronique comme variable a donné naissance à la théorie de la fonctionnelle de la 
densité (DFT), considérée comme l'approche la plus répandue et la plus polyvalente utilisée 


dans l'ère moderne de la matière condensée computationnelle. 
IL.3. Théorie de la fonctionnelle de la densité (DFT) 


Dans cette théorie, la densité électronique n(r) est utilisée pour résoudre l'équation de 
Schrüdinger qui est considérée comme la base de la DFT. En 1927, Thomas et Fermi [4,5] ont 
fait la première tentative d'utiliser la densité électronique pour obtenir plus d'informations sur 
les systèmes atomiques et moléculaires. Dans cette tentative, le modèle statistique quantique 
des électrons a été utilisé, car il ne considère que l’énergie cinétique des électrons. La 
contribution du noyau-électron et de l'électron-électron a été traitée classiquement. Dans ce 
modèle, ils ont déduit une expression simple pour obtenir l'énergie cinétique pour une densité 
de gaz d'électrons uniforme sans interaction, mais ont exclu l'échange et la corrélation des 


électrons. 
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En se basant sur un gaz d'électrons uniforme, Dirac [6] a étendu ce modèle en 
considérant le terme d'interaction d'échange et a modifié l'équation d'énergie cinétique. 
Cependant, les approximations ci-dessus (Thomas-Fermi et Dirac) ne décrivaient pas 


précisément l'électron dans le problème à plusieurs corps, ce qui a conduit à son échec. 
IL.3.1. Théorèmes de Hohenberg-Kohn 


La base solide de la DFT est apparue en 1964 sous la forme du formalisme de 


Hohenberg-Kohn [7], qui ont introduit leurs deux théorèmes suivants [8] : 


Théorème I : Pour tout système de particules en interaction dans un potentiel externe V,,4(T), 
ce dernier peut être déterminé de façon unique par la densité de particules de l'état 
fondamental n,(r) (en d'autres termes, le potentiel externe est une fonctionnelle unique de la 


densité). 


Théorème II : On peut définir une fonctionnelle universelle pour l'énergie E[n] en fonction de 
la densité n(r), valable pour tout potentiel externe V.,,4-(T), de sorte que l'énergie exacte de 
l'état fondamental du système est la valeur minimale globale de cette fonctionnelle et la 


densité n(r), qui minimise la fonctionnelle, est la densité exacte de l'état fondamental n,(r). 


En suivant ces deux théorèmes, l'énergie totale du système peut être écrite comme suit : 
EG) = FD + | Ver On dE (LS) 


alors que la fonctionnelle F[n(r)] a la forme suivante : 


F1] = TP] + JM] + Eraln®)|] (IL.6) 


Les termes T[n(r)], JIn(F)], et Esuln(r)] représentent l'énergie cinétique du système en 


interaction, le terme de Hartree (interaction classique entre les électrons) et la contribution 
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électrostatique (non-classique) provenant respectivement de l'auto-interaction (échange, c'est- 
à-dire la nature antisymétrique de l'électron) et des effets de corrélation des électrons, 


respectivement. 


Puisque le terme F[n(r)] ne dépend pas du V,,+(r), il reste le même pour tout système. 
Si F[n(r)] et n(r) étaient connus exactement, alors on peut facilement obtenir l'énergie de 
l'état fondamental et la densité dans un potentiel externe en minimisant la fonctionnelle, qui 
est une fonction de la densité tridimensionnelle. Cependant, les complications du système 
avec de nombreux électrons ne permettent pas d'obtenir la fonctionnelle F[n(r)] avec 
précision. Ainsi, les deux théorèmes cités ci-dessus ne fournissent pas avec précision la forme 
exacte de la fonctionnelle universelle. Heureusement, environ un an après l'article de 
Hohenberg-Kohn, Kohn et Sham [9] présentent une méthode simple pour effectuer des calculs 


DET, qui conserve la nature exacte de la DFT. 
IL.3.2. Formalisme de Kohn-Sham 


L'année suivante, en 1965, Kohn et Sham [9] ont introduit une approche pratique pour 
trouver la fonctionnelle inconnue en utilisant l’énergie cinétique d'un système sans interaction 
(Ts), construit à partir d'un ensemble d'orbitales, au lieu de l’énergie cinétique du système à 
plusieurs corps en interaction (T) et en gardant la même densité de l'état fondamental. Alors, 


pour un système sans interaction, le Tç est défini comme suit : 


1 OCC 
Ts = 3 D (WiV2p) AL.) 
i=1 
Kohn-Sham ont aussi défini la fonctionnelle d'énergie totale comme suit : 


Eln(?)] = Î Lund + Tr] + | BeS otre didf, + E[n(M] (IL.8) 
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où Vekr(r) et Ts sont les termes de potentiel externe et d'énergie cinétique, respectivement. Le 
troisième terme de l'équation (IL.8) est l'énergie électrostatique classique des électrons, connue 
sous le nom de terme de Hartree. E,. est le terme d'énergie d'échange-corrélation, et peut être 


défini comme suit : 
Exe = (TnP)] - Tsn)D + Enaln)] = Ten] + Evaln(”)] (IL.9) 


Ainsi, le terme E,. contient la véritable d’énergie cinétique, T+, et les contributions 
électrostatiques non classiques, E,4. La fonction d'énergie de Kohn-Sham définie dans 
l'équation (IL8) conduit à une équation de type Schrüdinger, connue sous le nom d'équation de 


Kohn-Sham. 


1 
Hs) = [592 +] = ei) (11.10) 


L'équation (11.10) montre que les corps non-interagissant se déplacent dans un potentiel 


effectif (VXs), et est défini comme suit : 
Vrs() = Vexe (7) + Va) + Vc() (IL 11) 


OÙ Voxt >» Vn et VLkc Sont respectivement le potentiel externe, le potentiel de Hartree et le 


potentiel d'échange-corrélation, respectivement. Ils sont donnés par les relations suivantes : 


nm) 
V, — ———<— d IL.12 
ré on 
5Exeln(F)] IL.13 
et Le — En ( ) 


42 


Chapitre II : Méthode de calcul 


S1 les fonctions propres et leurs valeurs propres correspondantes sont W; et £&;, la densité de 


l'électron de l'état fondamental peut être estimée comme suit : 
nG = D IE (11.14) 
i=1 


Maintenant, la densité électronique calculée dans l'équation (IL. 14) peut être utilisée pour 
estimer le nouveau potentiel effectif de manière cohérente. A partir des équations (11.7) et 


(IL.10), l’énergie cinétique d'un système sans interaction, peut être écrite comme suit ; 


OCC 


T[n(P)] = de - Î VA? (15) 


i=1 


et l'énergie totale peut être estimée en mettant l'équation (IL.15) dans l'équation (IL.8), soit : 


Em] = Ÿ a -2 Î 2 - Mr EMPAN + EMI (116 


Ir -#] 


Cette équation (IL.16) montre que la fonctionnelle d'énergie totale E[n(r)] est indépendante du 


potentiel externe V,,+(r). 
IL.4. Approximations de la fonctionnelle d'échange-corrélation 


Dans la section précédente, nous avons discuté de l'exactitude du formalisme de Kohn- 
Sham. Si le terme ÆE,. est connue, alors on peut calculer exactement l'énergie de l'état 
fondamental du système des particules en interaction. Bien que, la forme exacte de E,e est 
inconnue, on peut introduire une certaine approximation de E,-. Ainsi, la précision des 


calculs DET dépend uniquement de l'approximation choisie pour E,.. Différentes E,< ont été 
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introduites afin d'obtenir les énergies exactes de l'état fondamental et seront brièvement 


expliquées dans les sous-sections suivantes. 
IL.4.1. Approximation de la densité locale (LDA) 


La toute première forme de l'énergie de corrélation et d'échange a été proposée par 
Hohenberg et Kohn [7]. Dans cette proposition, ils ont considéré que la densité d'énergie 
d'échange-corrélation, £*{[n(#)], d'un système est identique à celui d’un gaz d'électrons 
uniforme de densité n(r). La fonctionnelle d'échange-corrélation pour LDA peut être définie 


comme suit : 


ELA [n(#)] = j ne (nd? (IL17) 


où el [n(r)] est la densité d'énergie d'échange-corrélation d'un gaz d'électrons uniforme avec 


une densité calculée localement, n(r), au point 7. C'est le mode le plus basique et 
fondamental de la fonctionnelle d'échange-corrélation, qui fonctionne de manière 
impressionnante pour de nombreux systèmes à l'état solide. €, est la contribution de deux 
termes ; le premier terme est l'échange (£,) et le deuxième est la corrélation (£.). Le terme €, 


peut être obtenu à partir de l'échange Hartree-Fock et dérivé par Dirac [6] : 


en] = (2) (U.18) 


Cependant, le mode précis du terme de corrélation est encore inconnu. Des valeurs très 
précises pour £&. ont été obtenues par Ceperley et Alder avec des simulations Monte-Carlo 
quantiques du gaz d'électrons homogène [10]. Sur la base des valeurs obtenues par Ceperley et 
Alder, un certain nombre d'expressions pour £. ont été proposées par différents auteurs 


utilisant des techniques plus avancées [11,12]. 
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La LDA semble bien fonctionner pour prédire les structures d'équilibre, les fréquences 
harmoniques et les propriétés magnétiques des systèmes moléculaires [13], cependant, elle 
sous-estime les paramètres du réseau [14] et a tendance à produire des énergies de liaison plus 


élevées (sur-liaison). 
IL.4.2. Approximation du gradient généralisé (GGA) 


Cette approximation a été introduite pour surmonter certaines limites de la LDA. Dans 
cette approximation, la densité électronique, n(r), et le gradient de la densité, V(F), ont été 
considérées pour approximer la densité d'échange-corrélation. L'énergie d'échange-corrélation 


de la GGA peut être estimée à partir de la densité d'énergie LDA, soit : 
LOTO) EE TOP TO) A TON ON" (119) 


ici F-[n(r), V(F)] est un terme supplémentaire, qui est connu comme la fonction 
d'amélioration. En 1992, Perdew et Wang ont approximé une forme d'énergie d'échange- 
corrélation, qui était libre du facteur d'amélioration [12,15]. Quelques années plus tard, cette 
forme a été simplifiée par Perdew, Burke et Ernzerhof (PBE). L’approximation GGA peut 
également être construite en utilisant le cadre de l'énergie de corrélation du gaz d'électrons et 
le terme dépendant du gradient [16]. La GGA a surmonté le problème de la sur-liaison de la 
LDA et a amélioré de nombreux aspects du matériau tels que les propriétés structurales, la 
transition de phase, la stabilité de phase et les énergies cohésives, etc. De plus, elle a 


considérablement amélioré les propriétés magnétiques itinérantes [13]. 


La LDA et la GGA décrivent avec succès certaines propriétés, comme mentionné ci- 
dessus, mais elles sous-estiment la bande interdite des semi-conducteurs et des isolants. Pour 
surmonter le problème des bandes interdites, d'autres approximations ont été introduites pour 


les estimer correctement. 
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IL.4.3. Fonctionnelles hybrides 


La fonctionnelle hybride [17] est une classe intéressante de fonctionnelles, qui inclut 


l'échange Hartree-Fock avec la GGA, dont la forme généralisée est : 


ER? = a(EËF — E$GA) + EGGA, (IL.20) 


où EXF est l'énergie d'échange de Hartree-Fock. Le coefficient «& indique la quantité de 
mélange d'échange exact, et ne peut être déduit à partir des méthodes de premiers principes. 


Donc, il est ajusté de manière semi-empirique. 


La méthode de Heyd-Scuseria-Ernzerhof (HSE) [17] a utilisé un potentiel de Coulomb 
blindé et une fonction d'erreur pour modifier la fonctionnelle d'échange-corrélation, afin 
d'estimer une partie de l'énergie d'échange, ce qui a permis d’améliorer l'efficacité des calculs, 


en particulier pour les systèmes métalliques. 


EP = ah (a) EE, Fo) +, 7 (a) HE FE (@), (11.21) 


où a, et © sont le mélange et le paramètre ajustable qui peut contrôler la gamme courte de 
l'interaction, respectivement. Les valeurs normalisées a = 0.25 et w = 0.20 se réfèrent à une 
méthode connue sous le nom de HSE06, qui s'est avérée donner des résultats raisonnables 
pour un certain nombre de systèmes. On peut obtenir la fonctionnelle hybride PBEO à partir 


de la fonctionnelle HSE avec « = 0. E. 


(w) représente la gamme courte de la 


fonctionnelle de Hartree-Fock, EX °*"* (w), Ex°""* 


(w), et EFPE(w) sont les composantes de 
gamme courte, longue de la fonctionnelle d'échange PBE, et les fonctionnelles de corrélation 


PBE, respectivement [16]. 
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IL.4.4. DET polarisée en spin 


Nous n'avons pas encore abordé la question de savoir comment la densité électronique 
dépend des électrons de spins différents et quel sera le formalisme DFT dans ce cas. Pour les 
systèmes avec des électrons non appariés, une extension de la DFT a été développée de la 
même manière que pour Hartree-Fock. Dans cette théorie, les densités d'électrons et les 
densités de spin sont les entités de base avec une densité de spin totale étant la différence entre 
les électrons de spin-up (spin majoriatire) et de spin-down (spin minoritaire) : p(r) = pr(r) — 
Pa(”). Par conséquent, la fonctionnelle d'échange-corrélation pour un système contenant des 


électrons de spins différents peut s'écrire comme suit : 


1 = 5 1 > ES > > 
|-5v + Vexe() + Î PC) Fate + VC, p)|DP ) = EPY? (F) (IL.22) 
2 


ici p(r) = pr(r) — pif”), et conduit le système à deux ensembles de fonction d'onde pour 


chaque état de spin. 
ILS. Solides périodiques 


Les calculs de structure électronique des systèmes solides sont confrontés à des 
problèmes de convergence en raison de la disponibilité d'un nombre infini d'électrons. Pour 
surmonter ce problème, on peut considérer la périodicité des solides. Dans un système à une 
seule particule, les électrons ont un potentiel, VX$, comme le montre l'équation de VX$ (voir 
équation (IL.10)), où Vxs(r + R), suit la périodicité du réseau. R représente la périodicité du 
réseau, similaire à celle du réseau de Bravais. Selon le théorème de Bloch [18], dans tout 
système périodique, la quantité ik est appelée moment cristallin du système, et le vecteur k 
est considéré comme un excellent nombre quantique qui administre la condition aux limites de 


la fonction d'onde (>), soit : 
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Da(F +R) = eRyP (es, 
où Ÿ> représente la fonction d'onde de Bloch : 
VO) = eue (7) (IL.24) 


où 4} (F) montre la fonction périodique du réseau (ux(r) = ux(r + R)). Afin de satisfaire le 
critère de Bloch de la condition aux limites, la fonction d'onde de la particule unique peut être 


étendue à un ensemble de base complet comme @;(r) : 


Var = > Cine Pix) (11.25) 


i 


où C;,x Sont les coefficients de l'expansion de Fourier. Maintenant, l'équation suivante peut 


être obtenue, en incluant l'équation (11.25) dans l'équation (IL.10) et en multipliant par (xl : 
à [loxlAslbx) — En lbxl6,x)] Cinr = 0 (11.26) 
j 


où les deux premiers termes représentent les éléments de la matrice hamiltonienne effective et 
de la matrice de chevauchement, respectivement. Les valeurs propres/coefficients d'expansion 


peuvent être obtenues en résolvant l'équation suivante : 


dett(é,:[#[ x) — Eux oxl PIC = (1.27) 
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IL.6. Ensembles de base : Ondes planes 


Un nombre considérable de difficultés numériques affectent encore la mise en œuvre de 
l'équation KS d'une seule particule en raison du comportement différent de la fonction d'onde 
dans différentes régions de l'espace, c'est-à-dire dans les régions du noyau et de la valence. 
Par conséquent, un ensemble de base complet est nécessaire pour décrire la fonction d'onde 


dans toutes les régions de l'espace. 


Il existe plusieurs choix possibles pour les ensembles de base en fonction du système 
étudié et de la précision requise comme : les ondes planes (PW), les ondes planes linéarisées 
augmentées (LAPW), les orbitales atomiques localisées (par exemple, LMTO : les orbitales 
linéaires Muffin-Tin), la combinaison linéaire des orbitales atomiques (LCAO), etc. Ici, nous 
allons mettre en lumière les ensembles de base des ondes planes, car la plupart des résultats de 


la thèse actuelle sont obtenus en utilisant des méthodes basées sur les ondes planes. 


La fonction périodique du réseau, 1 7 (r), peut être exprimée dans une série de Fourier 


de la forme suivante : 


hr = > Get” (1.28) 


_ 


G 


où c;g et G sont les coefficients d'expansion d'onde plane et le vecteur de réseau réciproque, 


respectivement, et G.r = 27m (m et Fr indiquent un entier et un vecteur de réseau cristallin, 
respectivement). À partir des résultats ci-dessus, la fonction d'onde de l'électron peut être 


développée sous la forme de combinaison linéaire d'ondes planes comme suit : 


A 1 Fer + 
Dr) = >. D ER O (II.29) 


= 
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où Ci et G sont les coefficient d'expansion de la fonction d'onde dans l'ensemble de base 


d'onde plane (el#+6)7) et les vecteurs réciproques du réseau, respectivement. Il est commode 


que les états soient normalisés et obéissent à une condition limite périodique dans un grand 


, — > : 1 
volume Q, que l'on laisse aller à l'infini. Par conséquent, le pré-facteur, a sert de facteur de 


normalisation, et le vecteur d'onde de Bloch est représenté par K. Par conséquent, l'équation 


KS dans la notation de l'état de Bloch s'écrit comme suit : 


h? 
(- V? +Vxs ©) Ÿ (F) = ErRŸ;r (F) (IL.30) 


2Me 


En incluant l'équation (11.29) dans l'équation (11.30), et en multipliant à gauche par eiE+6)F et 


en intégrant sur , nous obtenons l'équation aux valeurs propres de la matrice comme suit : 
ñ° + LA > A A E 
D Gr + 61833 + Vs (6 + Cox (ED) = exc (C) (11.31) 
G 


Dans la forme actuelle, l’énergie cinétique (premier terme) est diagonale, et Vxs est 
décrit sous la forme de ses transformés de Fourier. La solution de l'équation (1131) peut être 


2) avec les éléments donnés dans 


obtenue en diagonalisant la matrice hamiltonienne (H> 22 2 


les parenthèses du côté gauche. La taille de la matrice (somme sur G) est déterminée par une 


— 


: . . 2 2 2 . 
énergie de coupure appropriée, Ecut = re [k + Cal , et pourrait être plus grande pour les 
e 


systèmes ayant à la fois les électrons de valence et de noyau. On pense que c'est un problème 
grave, mais 1l pourrait être surmonté en utilisant l'approximation du pseudopotentiel, discuté 


dans la section suivante. 
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IL7. Les pseudopotentiels et les ondes augmentées projetées 


(PAW) 


Les fonctions d'onde atomiques sont orthogonales. La fonction d'onde de l'électron de 
valence doit osciller rapidement pour satisfaire la condition d'orthogonalité par rapport à la 
région du noyau. Ces oscillations sont confinées à la région centrale, qui n'est pas très 
importante pour l'interaction entre les atomes et qui nécessite de nombreuses ondes planes 
pour une description précise. Par conséquent, il pourrait être utile d'éliminer ces oscillations. 
Différents schémas ont été proposés pour traiter ce problème, deux des schémas les plus 


fréquemment utilisés sont décrits ci-dessous. 
IL.7.1. Les pseudopotentiels 


La méthode des pseudopotentiels est la méthode la plus utilisée. Elle comporte deux 
étapes : la première consiste à retirer les électrons du noyau et la seconde à remplacer les 
fonctions d'onde de valence par des pseudo-fonctions d’onde plus lisses. Les pseudo- 


fonctions d’ondes plus lisses nécessitent moins d'ondes planes pour leur description. 


Les électrons du noyau sont très compacts et étroitement liés à leurs noyaux hôtes et ne 
sont pas affectés lorsque les atomes interagissent entre eux. Par conséquent, il n'est pas grave 
de considérer les états du noyau comme fixes et de les éliminer des calculs. Le potentiel du 
noyau peut être combiné au potentiel nucléaire et ne traiter que les électrons de valence. Il 
faut garder à l'esprit que les fonctions d'onde associées aux électrons de valence (avec un 
moment angulaire différent) auront un potentiel différent en raison des conditions 
d'orthogonalité par rapport aux états du noyau. Par exemple, prenons l'exemple de l'état de 
valence 25 de l'atome de carbone qui a un potentiel différent de son état de valence 2p car les 
états 2s doivent être orthogonaux à l'état central 1s. Ainsi, le potentiel nouvellement obtenu 


est non-local (dans les coordonnées angulaires) et est connu sous le nom de pseudopotentiel. 


La deuxième étape consiste à remplacer le potentiel d'interaction du noyau ionique 
complet par un pseudopotentiel plus doux. Ainsi, les fonctions d'onde résultantes sont plus 


lisses et sans nœuds. De plus, les pseudopotentiels peuvent être encadrés en gardant leurs 
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propriétés de diffusion (déphasages) identiques au potentiel atomique réel dans l'énergie du 
premier ordre [19]. Cette dernière condition est satisfaite si la partie radiale de la fonction 
d'onde lisse (pseudo-fonction d'onde) est considérée comme identique à celle de la fonction 
d'onde réelle en dehors du rayon du noyau, 7, et si la dérivée de la pseudo-fonction d'onde est 
continue à rc, et à l'intérieur de la région du noyau, la norme de la fonction d'onde reste 
préservée simultanément. Un pseudopotentiel réel est construit dans l'autre sens, c'est-à-dire 
que l'on construit une pseudo-fonction d'onde qui satisfait ces conditions, et on inverse 


l'équation de Schrôüdinger pour obtenir le pseudopotentiel. 
IL.7.2. les ondes augmentées projetées (PAW) 


Une autre méthode souvent appliquée est la méthode de l'onde augmentée projetée 
(PAW) due au théorème de Bloch [20]. Par rapport à la méthode du pseudopotentiel, la 
méthode PAW préserve les électrons du noyau (dans un état gelé) et les fonctions d'onde de 
valence tout-électron. La fonction d'onde du noyau peut être prise à partir des estimations de 
tous les électrons des atomes considérés. La fonction d'onde lisse remplace la fonction d'onde 
réelle à l'intérieur du rayon du noyau (rc), ce qui est similaire à celui de la méthode du 
pseudopotentiel. La méthode PAW est basée sur une transformation linéaire de la pseudo- 


fonction d'onde Ÿ à la fonction d'onde tout-électron Ÿ suivante : 
W,)=16,)+ > Ib NGtn) - D 16) Gil) (1.32) 
i i 


avec Ÿ,, quantités variationnelles. L'indice i fait référence à une somme sur le site atomique, 
le moment angulaire et les énergies référencées. ® et ®, représentant l'ensemble des ondes 
tout-électron et pseudo-fonctions d’onde partielles qui correspondent parfaitement au rayon du 
noyau. Enfin, p; sont les fonctions de projecteur qui doivent être créées de manière à être 


doubles des fonctions d’onde partielles : 
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(Bild) = à; (IL.33) 


Ainsi, la méthode PAW, implémentée dans VASP, fonctionne directement avec la 
fonction d'onde de valence de tout-électron et les potentiels de valence tout-électron. Sa 


précision s'approche des méthodes tout-électron, et elle est plus exigeante en termes de calcul. 
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Chapitre III : 


Etude des propriétés structurales, 
électroniques et magnétiques de la 


monocouche Mol; dopée par des métaux 
de transition (TM : Sc, Ti, V, Cr et Mn) 


IIL.1. Introduction 


Dans ce chapitre, les propriétés structurales, électroniques et magnétiques de la 
monocouche Mol; pure et dopée par des métaux de transition (TM: Sc, Ti, V, Cr et Mn) ont 


été étudiées. La monocouche étudiée est dopée à une concentration de 12,5 %. 


III.2. Détails de calcul 


Nos calculs de premiers principes polarisés en spin ont été effectués à l'aide de la 
méthode des ondes augmentées par projecteur (PAW) [1], telle qu'implémentée dans le code 
Vienna ab- initio Simulation Package (V ASP) [2], basée sur la théorie de la fonctionnelle de la 
densité (DET). Pour tenir en compte des effets de corrélation, qui sont généralement 
significatifs dans les composés des métaux de transition, la méthode GGA + Ur introduite par 


Dudarev et al. [3] a été utilisée. Pour les atomes de Mo, nous avons choisi Ur = 3,0 eV, car 
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cela conduit à des résultats qui sont en bon accord avec la méthode HSEO06 [4], cependant, les 
valeurs de U; utilisées pour les atomes dopants des métaux de transition, elles ont été choisies 
selon des travaux antérieurs [5,7]. Dans tous les calculs, une énergie de coupure de 550 eV, 
pour l'expansion en onde plane des fonctions d'onde a été utilisée. Pour éviter les interactions 
entre deux images adjacentes, une supercellule 2x2x1 a été construite avec un vide de 15 À le 
long de la direction Z. Les constantes de réseau et les coordonnées atomiques ont été 
optimisées jusqu'à ce que la convergence des forces sur chaque atome soit inférieure à 0,01 
eV/À, tandis que le critère de convergence d'énergie était de 10% eV. Le schéma de 
Monkhorst-Pack centré en Gama [8] a été utilisé pour générer un maillage de points k 
permettant d'effectuer une intégration sur la première zone de Brillouin. Pour les calculs de la 
relaxation et de densités d'états, nous avons utilisé une grille de points k de 4x4xl et de 


15x15x1, respectivement. 
IIL.3. Propriétés structurales de la monocouche Mol; pure et dopée 


La monocouche Mol; est un membre de la famille des trihalogénures de métaux de 


transition, elle se présente sous une structure de couches en sandwich I-Mo-I, comme illustré 


Figure IIL.1 : (a) La partie supérieure (inférieure) représente la vue de dessus (vue de côté) de la 
structure cristalline de la monocouche Mol;. Une cellule unitaire de la monocouche Mol; se trouve à 
l'intérieur de la ligne solide en forme de losange. (b) Vue de dessus de la supercellule 2x2x1 de la 
monocouche Mol; dopée. 
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dans la figure IIL.1(a). Les ions magnétiques Mo** sont entourés dans la cellule unitaire par 
six ions I les plus proches voisins; chacun d'eux est lié à deux ions Mo** et organisé dans une 
coordination octaédrique, similaire à la monocouche Crl; synthétisée expérimentalement [9]. 
La monocouche Mol; dopée par les TM (TM = Sc, Ti, V, Cr et Mn) est obtenue à partir d’une 
supercellule 2xX2x1 qui contient 32 atomes (7 atomes de molybdène, 24 atomes d'iode et 1 


atome de TM), comme le montre la figure III. 1(b). 


D'après le tableau IIL.1, après une optimisation géométrique, la constante de réseau 
obtenue a pour la monocouche Mol; non dopée est de 7,241 À, ce qui est cohérent avec la 
valeur rapportée dans [4]. De plus, les constantes de réseau optimisées a, les longueurs de 
liaison d et les angles 0 des systèmes Mol; dopés par les TM sont comparables à ceux de la 
monocouche Mol; non dopée. Malgré les petites déformations trouvées, la structure cristalline 
des systèmes dopés est maintenue par rapport à la Mol; non dopée. Ensuite, nous avons étudié 
la stabilité relative des monocouches de Mol; non dopée et dopées par des TM en calculant 
l'énergie de formation £;,,. Pour la cellule unitaire de Mol; non dopée, nous avons obtenu 


une valeur de — 4,21 eV, calculée comme suit : 


Tableau IIL.1: Constante de réseau optimisée a, longueur de liaison TM-Mo d, angles de 
liaison TM-I-Mo & et énergie de formation £;,, dans les conditions riches en Mo et riches en I 
de la monocouche Mol; non dopée et dopée par les TM (TM = Sc, Ti, V, Cr et Mn). 


 . a (À) d (À) 0 (°) Eform (EV) 
Mo-riche 1-riche 
Pur 7,241 4,180 (dom) 94,269 (O-1-Mo) = s 
Sc-dopé 7,254 4,205 93,507 — 3,020 — 5,128 
Ti-dopé 7,242 4,186 94,33 — 1,498 — 3,605 
V-dopé 7,234 4,185 95,236 — 1,498 — 3,605 
Cr-dopé 7,227 4,146 94,72 — 0,463 — 2,570 
Mn-dopé 7,237 4,156 93,083 — 0,278 — 2,385 
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Ejorm = E(Molz) — E(Mo) — 3E(I2) (HIL.1) 


où E(Mol;) représente l'énergie totale de la monocouche Mol:, E(Mo) est l'énergie totale 
du Mo en état massif dans une cellule cubique centrée avec deux atomes de Mo, et E(L) 
représente l'énergie totale de la molécule IL. Cette grande valeur négative de E;,M indique la 
haute stabilité de la monocouche Mol; et correspond bien à la valeur de — 4,35 eV obtenue par 


la méthode HSEO06 [5]. 


Concernant les systèmes Mol; dopés par des métaux de transition (TM), les énergies de 


formation ont été calculées en utilisant la formule suivante [10] : 


EForm nn E(doped) — E(pristine) + n(Umo — Urm) (IIL.2) 


où les termes Eoped) €t Efpristiney ont référence aux énergies totales de la monocouche 
Mol avec et sans dopant, respectivement. Dans notre cas, n = 1 car un atome de Mo est 
remplacé par un atome de TM. h,,, et L,,, représentent les potentiels chimiques des atomes 
hôtes de Mo et des dopants TM, respectivement. ,,, est défini comme un ensemble de 
valeurs correspondant aux conditions de croissance riches en Mo ou riches en I. Ici, nous nous 
sommes uniquement concentrés sur deux conditions extrêmes, à savoir les conditions riches en 
Mo et les conditions riches en L. Pour une condition riche en Mo, h,,, est défini comme 


l'énergie d'un atome de Mo à l'intérieur de son réseau cubique centré stable, mais pour une 


_ E(Mols)-3E (2) 


condition riche en I, il est donné par la relation : wo : 


La valeur de u,,, de 


chaque dopant TM a été déterminée à partir de sa phase correspondante magnétique stable 
dans l’état massif. Les résultats présentés dans le tableau IT. 1 révèlent que tous les systèmes 
Mol dopés par des TM ont des énergies de formation négatives, ce qui indique que les 
processus de composition à partir des formes élémentaires sont des réactions exothermiques et 
que ces systèmes sont thermodynamiquement stables et ne se décomposent pas facilement une 
fois formés. Nos résultats montrent également que £;,, est plus faible dans les conditions 


riches en I, ce qui indique qu'il est plus énergétiquement efficace et plus facile d'incorporer des 
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atomes de TM dans Mol; dans des conditions expérimentales riches en I. De plus, le Mol; 
dopé par Sc dans les conditions riches en I est le plus stable, avec une énergie de formation 


minimale de — 5,128 eV. 


III.4. Propriétés électroniques et magnétiques de la monocouche 


Mol; pure 


Afin de voir l'effet du dopage sur la monocouche Mol, nous avons d'abord calculé ses 
propriétés électroniques et magnétiques à l'état pur. Les densités d'états (DOS) des spin-up et 
des spin-down représentées dans la figure IITL.2(b) sont asymétriques, indiquant ainsi la nature 
magnétique de la monocouche Mol; avec un moment magnétique total de 6,0 18 par maille 
unitaire (voir Tableau II2). De plus, nous avons confirmé la stabilité ferromagnétique de la 
monocouche Mol; en calculant la différence d'énergie magnétique entre les configurations 
ferromagnétique (FM) et antiferromagnétique (AFM), AË = Exy — Eaprm, et nous avons 
trouvé une valeur négative de  — 34,9 meV qui est en accord avec la valeur de — 36,5 meV 


obtenue par HSEO06 [5]. 


Tableau IIL2 : Moment magnétique total de la supercellule M,,, moment magnétique local 
des atomes de Molybdène M4, et d’Iode M;, gaps de spin-up Ep et gaps de spin-down 
E aps aps de charge E, et ES0C 


monocouche Mol; pure. 


ap Sans et avec couplage spin-orbite, respectivement, de la 


Système Mo: (ug) Muo(Ug) Mi(ls) Egap(eV) Egap(eV) Eg(eV) Ejap (eV) 


Pur 24,0 2,111 — 0,024 2,642 2,783 2,356 2,290 


De plus, nous remarquons d’après ces courbes que le maximum de la bande de valence 
(VBM) est principalement constitué des états de spin-up, tandis que le minimum de la bande 
de conduction (CBM) est principalement constitué des états de spin-down, ce qui fait de ce 
matériau un semi-conducteur ferromagnétique bipolaire (BFMS). En outre, les densités 
d’états partielles (PDOS) représentées dans la figure IIL.2(b) montrent que les états de VBM et 


CBM sont principalement constitués des orbitales Mo-4d et I-5p, ce qui indique l'existence 
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b) 


DOS (states/eV) 


Energy(eV) 


Figure IIL2 : (a) Structures de bandes polarisées en spin ; (b) Densités d’états totale et partielle 
(TDOS/PDOS) de la monocouche Mol; pure. 


d'une hybridation entre eux. D’autre part, la monocouche Mol; possède une symétrie 
ponctuelle C3, ce qui conduit à la division des orbitales Mo-4d en trois groupes : d,2, (dzy, 
dy2_y2), et (dx, dy,), comme le montre clairement la figure IIL.2(b). Aussi, les structures de 
bandes d’énergie représentées dans la figure IIL.2(a) révèlent que la monocouche Mol; à l'état 
initial possède une bande interdite directe de 2,356 eV située au point de symétrie M. Ces 


résultats sont en bon accord avec ceux obtenus par la méthode HSEO06 [41]. 


IIL.S. Propriétés électroniques de la monocouche Mol; dopée par 


des métaux de transition (TM : Sc, Ti, V, Cr et Mn) 


Les propriétés électroniques des systèmes Mol; dopés par des métaux de transition (TM) 
ont été étudiées en calculant les structures de bandes d’énergie polarisées en spin (Figure 
IIL3), les densités d'états totales (TDOS), les densités d'états partielles (PDOS) et l'effet SOC 
(couplage spin-orbite). Il ressort clairement des courbes PDOS illustrées aux figures IIL4, 
IILS et IIL.6 que les orbitales 3d de tous les dopants coexistent en un seul état a; (d,2) et deux 


états e1 (dxy/dy2_,2) et e2 (d;,/d,,) doublement dégénérés, car les atomes dopants sous le 
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champ cristallin présentent une symétrie C3. Dans les paragraphes suivants, nous discuterons 


plus en détail les résultats obtenus. 


Dans la monocouche Mol; dopée au scandium (Sc), un comportement semi-conducteur 
est obersvé dans les deux canaux de spin, avec un gap d’énergie de 2,267 eV le long du spin- 
up et 2,645 eV le long du spin-down situé au point de haute symétrie [', comme le montre la 
figure IIL.3(a). De même, un comportement semi-conducteur est également observé dans la 
monocouche Mol; dopée au titane (Ti), avec un gap d’énergie de 1,392 eV le long du spin-up 
et 2,4 eV le long du spin-down situé au point de symétrie l, comme le montre la figure 
IIL3(b). Cette diminution des valeurs de gap d’énergie dans les deux canaux de spin est 
probablement due à la présence de nouveaux états d'impuretés fortement localisés résultant 
initialement de l'hybridation entre les orbitales du dopant Ti et les atomes hôtes. En 
comparaison avec d'autres travaux, un comportement demi-métallique a également été observé 
en adsorbant des atomes de Sc ou de Ti sur une monocouche de Crl: [5], tandis que le dopage 
interstitiel de la monocouche VI à conduit à un comportement métallique ferromagnétique 
(FM) dans la monocouche dopée au Sc et à un comportement demi-métallique dans la 


monocouche dopée au Ti [11]. 


Le système Mol; dopé au vanadium (V) présente un comportement intéressant, avec une 
faible bande interdite d'énergie indirecte de 0,960 eV le long du canal de spin-up et une large 
bande interdite d'énergie directe de 2,497 eV le long du canal de spin-down, comme le montre 
la figure IIL3(c). En conséquence, la monocouche de Mol: dopée au V peut être classée 
comme un half semi-conducteur (HSC). Un comportement similaire a été observé dans une 
étude ultérieure sur le dopage interstitiel de VI: par des atomes de vanadium (V) [11]. Par 
contre, il a été démontré que l'adsorption d'atomes de vanadium (V) sur la monocouche Crl; 


entraîne un comportement demi-métallique [5]. 


Les structures de bandes d'énergie polarisées en spin de la monocouche de Mol; dopée 
au chrome (Cr), illustrées dans la figure IIL.3(d), révèlent l'existence d'une bande interdite 
d'énergie indirecte relativement étroite de 1,096 eV le long du canal de spin-up et d'une large 


bande d’énergie directe de 2,729 eV le long du canal de spin-down. Aïnsi, la monocouche 
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Figure IIL3 : Structures de bandes d’énergie polarisées en spin des monocouches Mol; dopées par les 
métaux de transition (TM : Sc, Ti, V, Cr et Mn), où les lignes rouges (bleues) représentent les canaux de 
Spin-up (spin-down). Le niveau de Fermi Er, réglé à zéro énergie, est indiqué par des lignes horizontales 

en pointillés noirs. 
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dopée au Cr présente le même comportement HSC que celle dopée au V, cependant, d'autres 
études ont montré que le dopage interstitiel de la monocouche VI; conduit à un comportement 
métallique [11], tandis que l'adsorption d'atomes de chrome (Cr) sur la monocouche Crl; 


entraîne un comportement demi-métallique [5]. 


On constate en particulier, que dans la monocouche Mol; dopée au Mn, les états 
d'impureté traversent le niveau de Fermi dans le canal de spin-up, alors qu'il existe un gap 
d'énergie appréciable de 2,619 eV dans le canal de spin-down, comme le montre la figure 
IIL3(e). Par conséquent, cette monocouche dopée présente un comportement demi-métallique 
et peut fournir un courant polarisé en spin de 100 % et peut ainsi être utilisée pour des 
applications de dispositifs de filtrage de spin. D'autre part, des études ont montré que le 
dopage interstitiel de la monocouche VI; avec du manganèse (Mn) présente des propriétés 
d’un semi-conducteur ferromagnétique bipolaire (BFMS) [11], tandis que l'adsorption 
d'atomes de manganèse (Mn) sur la monocouche Crl; conduit à des propriétés d’un semi- 
conducteur spin-gapless (SGS) qui se caractérise par un comportement semi-conducteur dans 
un canal de spin, tandis que dans le canal opposé, le VBM et le CBM se touchent au niveau de 


Fermi [5]. 


Pour donner plus de détails sur la contribution des orbitales des atomes aux gaps 
d'énergie créés et aux états d'impureté dans les structures dopées, nous avons tracé les courbes 
TDOS et PDOS pour chaque cas dopé. La figure IIL.4 révèle que le VBM est principalement 
due aux orbitales Mo-4d avec une faible contribution des états I-5p, tandis que dans le CBM, 
il y a des états d'impureté attribués principalement à l'hybridation entre les orbitales 3d de 
l'atome Sc et les orbitales 4d des atomes Mo les plus proches, avec une faible contribution des 
états 5p des atomes I les plus proches. Étant donné que le VBM et le CBM sont constituées 
par le même canal de spin, la monocouche Mol; dopée au Sc devient un semi-conducteur 
ferromagnétique plutôt que le BFMS qu'elle était dans son état pur. En revanche, dans le cas 
de la monocouche de Mol; dopée au Ti, on observe à partir des TDOS/PDOS dans la figure 
IIL.4 que la principale contribution au VBM et au CBM provient des orbitales Ti-3d, tandis 


qu'une faible contribution est due aux orbitales Mo-4d et I-Sp. 
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Figure IILA4 : Densités d'états totales et partielles (TDOS/PDOS) des monocouches Mol; dopées au 
Sc et au T1. 


67 


Chapitre III : Etude des propriétés structurales, électroniques et magnétiques de la 
monocouche Mol; dopée par des métaux de transition (TM : Sc, Ti, V, Cr et Mn) 


200 ' L H 7 Total - : + 
Ni Re | à Age 
2200 | | 1 

: 


T T T Mo-3d T 


DOS (states/eV) 
Non Non Non 
| 
< 
à 


= —. 
ec Re ee mans 


-2 L (9) 1 2 3 
Energy(eV) 
200 ke : ! H - | Total ÿ | : 
; ; Mo-3d —— 7 My! 1 
1 

A ñ l | ill IN 
| m— : 
o + Aeognshasr | 
-2 ! 1 2 1 
2 : T L 7 ] Cd —— F = | 

. _ 

> 0 

3 2 1 1 i 1 1 1 
5 2 7 L r ] Cry — T 2 
Z o _— sac | 
3 1 n si 


o 
Energy(eV) 


Figure IILS : Densités d'états totales et partielles (TDOS/PDOS) des monocouches Mol; dopées au 
V'et au Cr. 


68 


Chapitre III : Etude des propriétés structurales, électroniques et magnétiques de la 
monocouche Mol; dopée par des métaux de transition (TM : Sc, Ti, V, Cr et Mn) 


Concernant le cas du dopage au V, comme le montre la figure IILS, le VBM est 
principalement due aux orbitales Mo-4d et I-Sp, tandis que le CBM est principalement due aux 
orbitales V-3d. De plus, le gap d'énergie de spin-up a été réduit en raison de la formation de 
nouveaux états d'impureté hautement localisés dans la bande de conduction. Ces états 
d'impureté sont formés par le couplage des orbitales V-3d et Mo-4d ainsi que des orbitales 5p 
des atomes I, ce qui confirme l’existence d’une hybridation entre elles. Les PDOS du cas 
dopé au Cr dans la montrent que le VBM est formé par les orbitales Mo-4d et I-5Sp, tandis que 
le CBM provient de la contribution des orbitales Cr-3d et des orbitales I-Sp, ce qui indique 


également l'hybridation entre elles. 
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Figure IIL.6 : Densités d'états totales et partielles (TDOS/PDOS) de la monocouche Mol; dopée au 
Mn. 


Pour le cas dopé au Mn, représentée dans la figure IIL6, les PDOS montrent que les états 
d'impureté mentionné qui traversent le niveau de Fermi résultent de l'hybridation des orbitales 
Mo-4d, I-5p et Mn-3d. De plus, le canal de spin-up des deux orbitales e; (d,,/d,2_,2) et e2 


(dxz/dy,) doublement dégénéré est occupé, tandis que l'orbitale a; (d,2) est partiellement 
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occupé. Cela explique bien la valeur du moment magnétique local (4,544 u8) du Mn présentée 


dans le tableau IIL.3. 


Tableau IIL3 : Moment magnétique total de la supercellule M,,;, moment magnétique local 
du dopant Mrm ainsi que de ses atomes les plus proches voisins My et Mr, gaps de spin-up 
E ap et gaps de spin-down Ep; gaps de charge E, et ESC sans et avec couplage spin- 
orbite, respectivement, des monocouches Mol; dopées par TM (TM = Sc, Ti, V, Cr et Mn). 


Mru 


Système Mi (up) _. 
B 


Mio (us) Mile) Ejap(eV) Egap(eV) Eg(eV) ESF (eV) 


Sc-dopé 21,0 0,026 2,706  —0,020 2,267 2,645 2,207 2,005 
Ti-dopé 22,0 0,957 2,703 —0,029 1,392 2,400 1,392 2,256 
V-dopé 23,0 2,098 2,705  —0,028 0,960 2,497 0,960 1,587 
Cr-dopé 24,0 3,421 2,700 —0,078 1,096 2,729 1,096 2,005 
Mn-dopé 25,0 4,544 2,674  —0,034 0,000 2,619 0,000 0,000 


D'autre part, le couplage spin-orbite (SOC), considéré comme une forme du champ 
magnétique effectif "aperçu" par le spin de l'électron dans le référentiel du repos, est important 
dans la plupart des matériaux de la famille des trihalogénures de métaux de transition. Cela 
entraîne différents phénomènes, en particulier une réduction [5,12] et une ouverture 
(élargissement) des gaps d'énergie [13,14]. En plus, L’effet du SOC est étroitement liés au 
nombre d'électrons non-appariés et de la masse des atomes en question. Par conséquent, nous 
avons pris en compte le SOC pour étudier son effet sur les valeurs des gaps d’énergie, et 1l est 
possible de voir d’après le tableau IIL.3 qu'il y a une légère réduction des gaps d’énergie des 
monocouches Mol; non dopée et dopée au scandium et par conséquent, l'effet du SOC n'est 
pas très important dans ces deux systèmes. Cependant, dans les monocouches dopées au 
titane, au vanadium et au chrome, le SOC joue un rôle important dans l'élargissement des gaps 
d'énergie, tandis que dans le cas des monocouches dopées au manganèse, même avec l'effet du 


SOC, le système conserve son caractère demi-métallique. 
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IIL.6. Propriétés magnétiques de la monocouche Mol; dopée par 


des métaux de transition (TM : Sc, Ti, V, Cr et Mn) 


Il est évident que le dopage substitutionnel par des métaux de transition 3d peut 
influencer les propriétés magnétiques de la monocouche Mol. Les cinq systèmes dopés 
présentent un comportement magnétique, comme l'indique la différence entre les canaux de 
spin-up et spin-down des structures de bandes dans la figure IIL.3. Le moment magnétique 
total de la supercellule contenant un atome TM remplaçant un atome Mo est d'environ 21,0 LB, 
22,0 us, 23,0 8, 24,0 us et 25,0 up, respectivement pour les monocouches Mol; dopées au Sc, 
Ti, V, Cr et Mn, respectivement. Exceptionnellement, comme indiqué dans le tableau IIL3, le 
dopant Sc ne contribue à aucun moment magnétique, tandis que les dopants Ti, V, Cr et Mn 
ont des moments magnétiques locaux de 0,957 up, 2,098 up, 3,421 us et 4,544 up, 
respectivement. Parmi les cinq dopants TM, le dopant Mn a donc la plus grande contribution 
au moment magnétique total. De plus, les moments magnétiques négatifs des atomes d'iode 
() voisins indiquent l'existence d'un couplage antiferromagnétique entre eux, ainsi qu'entre les 
dopants et les atomes de molybdène (Mo). Il est également à noter que la valeur négative la 


plus élevée se trouve dans la monocouche Mol; dopée au chrome (Cr). 


Pour obtenir une vue détaillée sur la distribution des moments magnétiques des systèmes 
Mol dopés aux TM, nous avons tracé les isosurfaces de densités de spin telles qu'illustrées 
dans la figure IIL.7. Dans tous les systèmes dopés, la densité de spin est dominée par les 
électrons de spin "up" et elle se trouve principalement localisée autour des atomes de Mo et 
des dopants TM, à l'exception du dopant Sc qui ne présente aucune contribution à la densité de 
spin. De plus, on remarque que les atomes d'iode voisins sont couplés de manière 
antiferromagnétique aux atomes de Mo et aux dopants TM et ce couplage est plus important 
dans le cas du dopant Cr. En particulier, les dopants Cr et Mn ont une distribution de densité 
de spin "up" quasi-sphérique autour d'eux, et ils induisent une densité de spin plus élevée sur 
les atomes voisins par rapport aux autres dopants. Ces résultats sont en accord avec ce qui a 


été discuté précédemment. 
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Figure IIL7 : Densités de spin des monocouches Mol; dopées au a) Sc, b) Ti, c) V, d) Cret 
e) Mn, où les isosurfaces jaunes (cyan) représentent le spin-up (spin-down ) densité (+0,0022 
23 


Comme le robuste ferromagnétisme et la haute température de Curie sont essentiels pour 
les applications spintroniques du monde réel, nous avons calculé la différence d'énergie 
magnétique (AE) entre les configurations FM et AFM, ainsi que la température de Curie dans 
l’approximation de champ moyen (TF4). Nous pouvons observer à partir des valeurs 
résumées dans le tableau IIL4 qu’à l'exception du cas dopé au scandium, tous les dopants 
présentent des valeurs négatives de AE plus élevées que le cas d'origine, indiquant un 
renforcement de la stabilité ferromagnétique. Bien que l'approximation du champ moyen a 
tendance à surestimer la température de Curie, elle peut donner une estimation de 


l'augmentation de la température de Curie due au dopage et fournir ainsi une base théorique 
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MFA 
Te 


pour de futures études expérimentales. Nous avons donc estimé en utilisant la relation 


2 . 
suivante [15-18] : k&TAMFA = AE, où k4 est la constante de Boltzmann. On peut constater 


que les valeurs rapportées dans le tableau IIL4 sont supérieures à celles obtenues pour la 
monocouche Mol; d'origine et bien supérieures à celles obtenues dans certains matériaux 
bidimensionnels purs et dopés avec des métaux de transition tels que le MoS; [19], le silicane 
[16], le SnS [20] et le Crl; [5]. De manière encore plus intéressante, les valeurs élevées de 
TMFA observées dans les cas dopés au titane et au manganèse, dépassant 600 K, prédisant ainsi 
qu'ils pourraient avoir une température de Curie élevée même dans la réalité, ce qui en fait des 
candidats très prometteurs pour les applications de dispositifs bidimensionnels de stockage de 


données. 


Tableau IIL.4 : La différence d'énergie magnétique entre les configurations FM et AFM, AE, 
ainsi que la température de Curie dans l’approximation du champ moyen TAFA des 
monocouches Mol; pure et dopées par des métaux de transition (TM = Sc, Ti, V, Cr et Mn). 


Système AE (eV) DFA) 


Pur — 0,03498 270,617 
Sc-dopé — 0,03465 268,064 
Ti-dopé — 0,07818 604,827 
V-dopé — 0,03773 291,892 
Cr-dopé — 0,04145 320,671 
Mn-dopé — 0,08336 644,901 


IIL.7. Conclusion 


En résumé, nous avons utilisé des calculs de premiers principes basés sur l'approche 
GGA+U:r pour étudier les propriétés structurales, électroniques et magnétiques des 
monocouches Mol, à la fois pure et dopées avec des métaux de transition (TM = Sc, Ti, V, Cr 
et Mn). Nos résultats montrent que le dopage substitutionnel avec des métaux de transition 3d 


est une technique efficace pour améliorer les propriétés spintroniques de la monocouche Mol. 
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Les monocouches dopées avec Sc et Ti présentent des propriétés semi-conductrices 
ferromagnétiques (FM), tandis que les monocouches dopées avec V et Cr présentent des 
propriétés half semi-conductrices. Ces quatre systèmes possèdent différents gaps d’énergie 
qui sont sujettes à l'effet du couplage spin-orbite (SOC), à l'exception de la monocouche dopée 
avec Mn, qui présente un comportement demi-métallique avec une polarisation de spin 


complète et un magnétisme renforcé. 


De plus, la stabilité ferromagnétique et la température de Curie dans l’approximation du 


TMFA ont été considérablement améliorées grâce au dopage substitutionnel, en 


champ moyen 
particulier dans le cas de la monocouche dopée avec Ti et Mn. Par conséquent, les 
monocouches dopées avec Sc, Ti, V, Cr et Mn présentent un énorme potentiel pour les 
applications spintroniques de semi-conducteurs de nouvelle génération avec différentes 


propriétés électroniques et magnétiques. 
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Chapitre IV : 


Etude des propriétés structurales, 
électroniques et magnétiques de la 


monocouche VCT,; dopée par des métaux 
de transition (TM : Sc, Ti, Cr, Mn et Fe) 


IV.1. Introduction 


Dans ce chapitre, nous avons étudié les propriétés structurales, électroniques et 
magnétiques de la monocouche VC, à la fois pure et dopée par des métaux de transition (TM 


: Sc, Ti, Cr, Mn et Fe) avec la même concentration indiquée précédemment. 


IV.2. Détails de calcul 


Comme indiqué dans le chapitre précèdent, nous avons utilisé l'approche des ondes 
augmentées par projecteur (PAW) [1], implémentée dans le code Vienna Ab-initio Simulation 
Package (VASP) [2] et les calculs de premiers principes polarisés en spin ont été réalisés dans 
le cadre de la théorie de la fonctionnelle de la densité (DFT). Pour décrire les effets de 
corrélation, qui sont souvent considérables dans les systèmes contenant des métaux de 


transition 3d, nous avons appliqué la correction Hubbard "U" dans la méthode GGA + Ur 
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invariante par rotation, développée par Dudarrev et al. [3]. Dans cette approche, la valeur 
Ueg = U — J est d'une importance particulière par rapport aux paramètres individuels U et J. 
Concernant les atomes V, nous avons adopté la valeur U,; = 3,35 eV, déterminée de manière 
auto-cohérente par [4]. Quant aux valeurs de U,; utilisées pour les dopants, elles ont été 
sélectionnées en accord avec des publications antérieures [5-7]. De plus, nous avons employé 
une énergie de coupure de 550 eV pour la base d'ondes planes dans tous les calculs. Pour 
éviter les interactions entre les couches adjacentes, nous avons construit une supercellule de 
dimensions 2x2x1 avec un vide de 15 À le long de la direction Z. Les constantes de réseau et 
les positions atomiques ont été optimisées jusqu'à ce que les forces de Hellmann-Feynman sur 
chaque atome soient inférieures à 0,01 eV/À, tandis que le critère de convergence énergétique 
a été fixé à 10° eV. Pour l'optimisation structurale et les calculs électroniques, nous avons 
utilisé respectivement une grille de points k de 4x4x1 et de 15x15x1, centrée sur Gamma 


selon la méthode de Monkhorst-Pack [8]. 
IV.3. Propriétés structurales de la monocouche VCTI; pure et dopée 


La monocouche VC appartient à la famille des trihalogénures de métaux de transition. 
Elle présente une structure en couches similaire à celle du Mol, comme illustré dans la figure 


IV.1. Nous avons appliqué la même procédure d’étude précédente du Mol; pour l’étude des 


Figure IV.1 : Vues de dessus et de côté de la monocouche VCI; (a) pure et (b) dopée par des 
éléments de transition (TM). 
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propriétés structurales de la monocouche VC et les résultats obtenus sont résumés dans le 
tableau IV.1. La valeur trouvée de la constante de réseau optimisé est de 6,325 À, ce qui est 
en bon accord avec les résultats antérieurs de 6,28 À [4] et 6,25 À [9]. De plus, nous 
remarquons d’après les valeurs rapportées dans ce tableau, que l'énergie de formation Ep est 
négative et a une valeur de — 12,27 eV par maille unitaire et de — 1,53 eV par atome, 
confirmant ainsi la stabilité thermodynamique de la monocouche VCI:, tel que trouvé par la 


méthode HSE06 (-— 1,37 eV/atom) [10]. 


Tableau IV.1 : Constante de réseau optimisée, a, longueur de liaison TM-V, d, angles de 
liaison TM-CI-V, 0 et énergie de formation E;, dans des conditions riches en V et riches en 
CI des monocouches VC pure et dopées aux éléments de transition (TM = Sc, Ti, Cr, Mn et 
Fe). 


ee a (À) d (À) 0 (°) Eform (eV) 
V-riche Cl-riche 
Pur 6,325 3,651 (dy.y) 98,244 (0y.c1.v) = 
Sc-dopé 6,301 3,698 96,68 — 5,049 — 11,186 
Ti-dopé 6,281 3,610 96,54 — 3,198 — 9,335 
Cr-dopé 6,245 3,554 95,37 1,702 — 7,839 
Mn-dopé 6,289 3,579 93,56 — 0,834 6,971 
Fe-dopé 6,257 3,628 97,18 — 0,698 — 6,835 


Pour le dopage substitutionnel de la monocouche VC, cinq éléments de transition (Sc, 
Ti, Cr, Mn et Fe) ont été choisis. Toutes les énergies de formation des systèmes dopés aux 
éléments de transition sont négatives, ce qui confirme leur stabilité thermodynamique. Il est à 
noter que la plus faible énergie de formation est observée dans des conditions riches en chlore 
que celles riches en vanadium, suggérant que l'insertion des éléments de transition dans la 
monocouche VC: est énergétiquement plus favorable dans de telles conditions 
expérimentales. Plus encore, Le dopage avec le Sc se révèle être le plus stable sous ces 


conditions, avec une énergie de formation minimale de — 11,186 eV. 
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IV.4. Propriétés électroniques et magnétique de la monocouche 


VC pure 


Dans un premier temps, afin de comprendre l'influence du dopage substitutionnel sur la 
monocouche VC, nous avons calculé les propriétés électroniques et magnétiques de celle-ci 
dans son état pur. La différence observée dans la distribution de la densité d'états (DOS) entre 
les deux canaux de spin, illustrée dans la figure IV.2(a), indique que la monocouche VCI est 
magnétique avec un moment magnétique total de 4,0 u8 par maille unitaire (voir tableau IV.2). 
Il en découle que la contribution prédominante au moment magnétique total provient des deux 
atomes de vanadium, avec 2 8 par chacun, tandis que la contribution des atomes de chlore 
environnants est négligeable. De plus, nous avons calculé la différence d'énergie magnétique 
ÂËE et l'avons trouvée égale à — 34,6 meV, confirmant la stabilité de la monocouche VCI; dans 
son état ferromagnétique, comme indiqué par He et al. [4]. Selon les structures de bandes 
polarisées en spin, illustrées dans la figure IV.2(b), la monocouche VC]; pure présente un gap 


d’énergie nul le long du canal de spin-up, et un large gap d’énergie indirect de 4,635 eV le 


(a) 


Total — 
40 L VE0d J 
CI-3p — 
20 | 
-40 | | 
T T 


Energy(eV) 
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DOS (states/eV) 
Energy (eV) 


Figure IV.2 : (a) TDOS et PDOS ; (b) Structures de bandes de la monocouche VC: pure. Le niveau de 
Fermi Er, placé à zéro énergie, est indiqué par des lignes en pointillés noirs. 
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long du canal de spin-down. De plus, l'existence d'un cône de Dirac dans le canal de spin-up 
permet de classer la monocouche VCH d'origine comme un demi-métal de Dirac, comme a été 


suggéré par [4]. 


Tableau IV.2 : Moment magnétique total de la supercellule M, moment magnétique local 
des atomes de Vanadium M, et de Chlore et Mc, gaps de spin-up E Dan et gaps de spin-down 


E ins respectivement, de la monocouche VCL pure. 
| Mio (us) My(usg) Ma(us) Egap (@V) 
Système : l 
Egap (CV) Egap (CV) 
Pur 16,0 2,060 — 0,072 0 4,635 


En général, étant donné que la monocouche VC est placée dans un champ cristallin 
octaédrique, les orbitales V-3d se sépareraient en un triplé inférieur t2, (dx,, d,, et dx) et un 
doublet supérieur e, (d,2 et d,2_,2). Selon Yang et al. [11], dans la symétrie P31M, le triplet 
t2g Se diviserait ensuite en un singlet a et un doublet €, et en raison de la rotation du cristal, le 
sous-état a le plus bas est d,2 et le doublet e, est d,., et d,2_,2. Cela correspond bien aux 
DOS/PDOS représentés dans la figure IV.2(a) et il est également clair que les orbitales 3d 
dans le plan (d,2_,2 et d,,) contribuent principalement aux états autour du niveau de Fermi 
dans le canal de spin-up (Figure IV.2(a)). De plus, le maximum de la bande de valence 
(VBM) est constitué par les états Cl-3p, tandis que le minimum de la bande de conduction 
(CBM) est formé par les états V-3d le long du canal de spin-down. Ces résultats concordent 


assez bien avec des études antérieures [4, 12, 13]. 


IV.5. Propriétés électroniques de la monocouche VCI; dopée par 


des métaux de transition (TM : Sc, Ti, Cr, Mn et Fe) 


Les propriétés électroniques de tous les systèmes VCI; dopés par des métaux de 
transition (TM) ont été explorées en procédant au calcul des structures de bandes polarisées en 
spin (Figure IV.3), ainsi qu'en évaluant les TDOS et PDOS. Veuillez noter que les PDOS des 


dopants sont tracées avec des échelles différentes afin de révéler les détails de chaque DOS 
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projetée. Il est possible de voir, à partir de la figure IV.3, que les propriétés électroniques de 
la monocouche VCI3 dopée par des métaux de transition dépendent fortement du type du 
dopant introduit. Cela est certainement en étroite relation avec la différence entre les électrons 
de valence du dopant et de l'atome de vanadium substitué comme l'augmentation ou la 
diminution du nombre d'électrons de valence qui s'accompagne d'un déplacement du niveau de 
Fermi et de la création de bandes d'impuretés provenant de l'hybridation entre les états des 
atomes hôtes et du dopant. Il est également possible de remarquer, à partir de la figure IV.3, 
que l'effet du dopant Fe est significatif sur les deux canaux de spin par rapport au reste des 
dopants, où l'effet sur le canal de spin-down est très faible. Cela peut s'expliquer par la grande 


différence de nombre d'électrons de valence entre le dopant Fe et les atomes V [14]. 


Comme indiqué dans les trois premiers graphes de structure de bandes polarisées en 
spin, le dopage de la monocouche VC avec les atomes Sc, Ti ou Cr entraîne la création de 
bandes interdites avec différents gaps d’énergie dans le canal de spin-up et donc un 
changement de la structure électronique de la monocouche VCI; du demi-métal de Dirac, qui 
était dans le cas pur, à un semi-conducteur ferromagnétique à l'état dopé. On peut également 
voir sur la figure IV.3(a) et la figure IV.3(b) que les cas dopés avec Sc et Ti sont des semi- 
conducteurs avec des gaps d’énergie indirects le long du canal de spin-up et sont égaux à 
1,863 eV et à 0,811 eV, respectivement. En ce qui concerne le canal de spin-down, la 
monocouche VCI; dopée avec Sc possède un gap d’énergie indirect égal à 3,602 eV, tandis 
que dans le cas dopé Ti, il existe une bande interdite directe de 3,568 eV située au point de 


haute symétrie M. 


D'autre part, la monocouche VCI3 dopée au Cr présente des bandes interdites directes 


situées au point de haute symétrie M pour les deux canaux de spin, où la valeur de la bande 
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Figure I1V.3 : Structures de bandes d’énergie polarisées en spin des monocouches VCI; dopées par les 
métaux de transition (TM = Sc, Ti, Cr, Mn et Fe), où les lignes rouges (bleues) représentent les canaux 
de spin-up (spin-down). Le niveau de Fermi Er, réglé à zéro énergie, est indiqué par des lignes 
horizontales en pointillés noirs. 


85 


Chapitre IV : Etude des propriétés structurales, électroniques et magnétiques de la 
monocouche VCI; dopée par des métaux de transition (TM : Sc, T1, Cr, Mn et Fe) 


interdite du canal de spin-up est égale à 1,721 eV et celle du canal de spin-down est égale à 


3,540 eV, comme illustré dans la figure IV.3(c). 


Pour connaitre la contribution des orbitales atomiques aux gaps d’énergie dans les 
systèmes monocouches VCI3 dopés aux métaux de transition (TM), nous avons tracé les 
figures PDOS (Figures IV.4-IV.6) pour chaque cas de dopage. Nous constatons, d’après ces 
figures, que la dégénérescence des orbitales 3d, qui était présente dans l'état pur de la 
monocouche VCB, est levée après le dopage. Ceci est probablement attribué à la rupture de 
symétrie de la structure causée par la présence du dopant, qui a créé des dislocations et des 
déformations dans le réseau (Tableau IV.1). Un comportement similaire a été rapporté dans 
des études antérieures où la présence de dopants a provoqué une levée de dégénérescence des 
orbitales d [15, 10]. Comme le montre la figure IV.4, les VBM et CBM sont principalement 
formés par les états V-3d dans les deux cas de monocouches VCI3 dopées au Sc et au Cr. En 
revanche, dans le cas dopé au Ti, la plus grande contribution au VBM provient des états Ti-3d 
(d 


xy» dyz) du dopant, tandis que le CBM est principalement contribué par les états V-3d, 


comme illustré dans la figure IV.5. 


D'autre part, les TDOS/PDOS du cas dopé au Mn, représentées dans la figure IV.S, 
révèlent que les VBM et CBM sont dominés par les états Cl-3p et Mn-3d (d,2_,2, dxz), 
impliquant une forte hybridation p-d autour du niveau de Fermi. En plus de la fourniture d’un 
courant polarisé en spin, la nature semiconductrice spin-gapless (SGS) de la monocouche 
dopée au Mn peut être ajustée par le champ magnétique et la contrainte pour diverses 


applications en spintronique, optique et capteurs [16]. 


En ce qui concerne la structure semiconductrice ferromagnétique bipolaire (BFMS) du 
cas dopé au Fe, les TDOS/PDOS présentées dans la figure IV.6, montrent que le VBM est 
constitué uniquement par les états V-3d, tandis que le CBM est dominé par les états Fe-3d du 
dopant et les états Cl-3p. Ainsi, le couplage des états Fe-3d et des états 3p des atomes de 
chlore environnants conduit à la formation de nouveaux états d'impuretés hautement localisés 
dans la bande de conduction, indiquant une hybridation p-d significative entre eux. La nature 


BFMS de VC; dopé au Fe offre la possibilité de contrôler électriquement l'orientation du spin 
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Figure IV.4 : Densités d'états totales et partielles (TDOS/PDOS) des monocouches VCl: dopées au 
Sc et Cr. 
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Figure IV.5 : Densités d'états totales et partielles (TDOS/PDOS) des monocouches VCI: dopées au 
Ti et Mn. 
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Figure IV.6 : Densités d'états totales et partielles (TDOS/PDOS) de la monocouche 
VCH dopée au Fe. 


des porteurs par application de tension, ce qui permet à cette monocouche dopée d'être utilisée 
pour la fabrication de plusieurs dispositifs de spintronique tels que des vannes de spin à effet 


de champ et des filtres de spin à effet de champ bipolaires [17]. 


IV.6. Propriétés magnétiques de la monocouche VCI; dopée par 


des métaux de transition (TM : Sc, Ti, Cr, Mn et Fe) 


Nous avons étudié l'influence des dopants TM sur les propriétés magnétiques de la 
monocouche VC et les résultats obtenus sont résumés dans le tableau IV.3 précédent pour 
l'état pur et chaque état dopé. Comme le montre les valeurs rapportées dans tableau IV.3, 
lorsqu'un atome de V est substitué par un atome de TM, le moment magnétique total de la 
supercellule devient égal à 14,0 u8, 15,0 u8, 17,0 uB, 18,0 ug et 19,0 u8 pour les monocouches 


dopées au Sc, Ti, Cr, Mn et Fe, respectivement. Le dopant Sc ne contribue à aucun moment 
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magnétique local, tandis que le dopant Ti contribue à un moment magnétique local de 0,892 


LB, réduisant ainsi le moment magnétique total de la supercellule qui était de 16,0 us à l'état 


Tableau IV.3 : Moment magnétique total de la supercellule M,,;, moment magnétique local du 
dopant Mry ainsi que de ses atomes voisins les plus proches My et Mc, gaps de spin-up E de 
et gaps de spin-down E D respectivement, des monocouches VCI; dopées par (TM = Sc, Ti, 


Cr, Mn et Fe). 


| Mi (Us) Mru(ue) My(ug) Ma(u) Egap (CV) 
Système 
Evan (CV) Esp (CV) 

Sc-dopé 14,0 0,004 1,941 — 0,013 1,863 3,606 
Ti-dopé 15,0 0,892 1,935 — 0,022 0,811 3,568 
Cr-dopé 17,0 3,068 1,938 — 0,043 1,721 3,540 
Mn-dopé 18,0 4,495 1,937 — 0,033 0,009 3,401 
Fe-dopé 19,0 4,106 1,938 0,095 1,624 2,343 


pur. Cependant, les dopants restants entraînent une augmentation du moment magnétique de 
la supercellule. Le dopant Mn contribue le plus avec un moment magnétique local de 4,495 
us, tandis que le plus grand moment magnétique total est trouvé dans la monocouche VCI; 
dopée au Fe. En particulier, la valeur positive du moment magnétique local de l'atome de CI 
dans le cas dopé au Fe indique que les atomes de CI voisins sont couplés ferromagnétiquement 
avec le dopant Fe, tandis que le couplage est antiferromagnétique dans les autres cas de 


dopage, comme le montre les moments magnétiques locaux négatifs des atomes de CI. 


Afin d'obtenir une vue détaillée de la distribution des moments magnétiques dans les 
systèmes VCI; dopés en TM, les isosurfaces de densité de spin sont représentées dans la 
figure IV.7, et il est possible de constater que ce sont les électrons de spin-up qui dominent la 
densité de spin de tous les systèmes dopés. De plus, on constate que la densité de spin est 
principalement localisée autour des atomes de V et des dopants TM, à l'exception du dopant 
Sc qui ne contribue à aucune densité de spin. De manière exceptionnelle, les atomes de CI 
voisins sont découverts comme étant couplés ferromagnétiquement au dopant Fe, renforçant 


ainsi le ferromagnétisme de la monocouche VC; et augmentant son moment magnétique total. 
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En revanche, les atomes de CI voisins dans les autres systèmes dopés sont couplés 
antiferromagnétiquement aux atomes de V et aux autres dopants TM restants. Ces résultats 


sont cohérents avec ce qui a été discuté précédemment. 


Figure IV.7 : Densités de spin des monocouches VCI; dopées aux métaux de transition, où 
les isosurfaces jaunes (cyan) représentent le spin-up (spin-down ) densité (+0,003 e/A3). 


IV.7. Conclusion 


En conclusion, nous avons calculé et discuté les propriétés structurales, électroniques et 
magnétiques des systèmes monocouches VCI; pur et dopés avec des TM (TM = Sc, Ti, Cr, Mn 
et Fe). Nous avons montré clairement la possibilité d’améliorer les propriétés spintroniques de 


la monocouche VCI, en la dopant avec des métaux de transition 3d. 
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Les monocouches dopées au Sc et au Ti révèlent un comportement semiconducteur 
ferromagnétique avec des gaps d’énergie indirects dans le canal de spin-up, tandis que la 
monocouche dopée au Cr présente des gaps d’énergie directs dans les deux canaux de spin, 
conduisant à une monocouche semiconductrice ferromagnétique à gap directe. Plus 
intéressant encore, la monocouche VCIl; dopée au Mn entraîne un comportement 
semiconducteur spin-gapless (SGS), tandis que celle dopé au Fe s'avère être un 


semiconducteur ferromagnétique bipolaire (BFMS). 


De plus, d'importants moments magnétiques de 4,25 18, 4,5 ue et 4,75 uB par maille 
unitaire sont observés respectivement dans les cas dopés au Cr, Mn et Fe, et de manière 
exceptionnelle, un couplage ferromagnétique est observé entre le dopant Fe et les atomes de 


CI voisins, renforçant la stabilité ferromagnétique de la monocouche VC. 


Ces résultats ouvrent une nouvelle voie pour améliorer les propriétés spintroniques des 
matériaux 2D à magnétisme intrinsèque. Notre étude peut être étendue à l'utilisation d'autres 
éléments TM 4d et 5d, et ouvre la voie au développement de dispositifs électroniques et 


spintroniques à haute performance. 
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Conclusion générale 


Dans le contexte de cette thèse, nous avons mené une étude de l’effet du dopage 
substitutionnel sur les propriétés structurales, électroniques et magnétiques des monocouches 
Mol et VC qui sont des membres de la famille des trihalogénures de métaux de transition MX3. 
Pour atteindre notre objectif, nous avons utilisé une approche de calcul ab-initio basée sur la 
théorie de la fonctionnelle de la densité (DFT), en employant la méthode des ondes augmentées 


par projecteur (PAW) et l'approche GGA+Urr telle qu'implémentée dans le code VASP. 


Nous avons commencé par l’étude des propriétés structurales qui définissent l'état 
fondamental des systèmes étudiés, notamment le paramètre de réseau calculé à l'équilibre, qui 
correspond bien aux données expérimentales et théoriques disponibles. Ensuite, nous avons 
examiné la stabilité thermodynamique des monocouches Mol; pure et dopées en calculant leurs 
énergies de formation sous des conditions riches en Mo et riches en I. Les résultats obtenus ont 
révélé que l'introduction des métaux de transition (TM) dans les monocouches Mol: conduisait à 
des énergies de formation négatives, signifiant ainsi des réactions exothermiques et démontrant 
une stabilité thermodynamique de ces systèmes après leur formation, de telle sorte que 
l'incorporation d'atomes de TM dans le Mol est plus efficace et favorable sur le plan énergétique 
dans des conditions expérimentales riches en I que celles riches en Mo. Parmi les différents 
systèmes dopés, le Mol; dopé au scandium a présenté la plus grande stabilité, avec la plus faible 


énergie de formation indiquant son potentiel d'adéquation pour diverses applications. 


Concernant les propriétés électroniques et magnétiques de la monocouche Mol, le dopage 
substitutionnel a été fait par les métaux de transition 3d (TM = Sc, Ti, V, Cr et Mn), avec une 
concentration de 12,5%. L’analyse des figures TDOS/PDOS a révélé la création de nouveaux 
états d'impureté autour du niveau de Fermi responsables de la modification des propriétés 
électroniques et, par conséquent, les gaps d’énergie des systèmes dopés avec des métaux de 
transition. Ces états d'impureté sont attribués à l'hybridation p-d entre les dopants et les atomes 
hôtes (Mo et I). Nos résultats montrent clairement que lorsque les monocouches sont dopées 
avec du Sc et du T1, elles présentent des propriétés semi-conductrices ferromagnétiques (FM), 
tandis que celles dopées avec du V et du chrome Cr, elles présentent des propriétés half semi- 
conductrices. Ces quatre systèmes possèdent des gaps d'énergie différents qui sont influencés par 
le couplage spin-orbite (SOC), à l'exception de la monocouche dopée au Mn, qui se comporte de 


manière demi-métallique avec une polarisation de spin complète et un magnétisme renforcé. De 
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plus, il est à noter que le dopage substitutionnel améliore considérablement la stabilité 


ferromagnétique et la température de Curie dans l'approximation du champ moyen TMFA, 


en 
particulier dans le cas de la monocouche dopée avec du Ti et du Mn. Cette amélioration de la 
stabilité ferromagnétique se manifeste notamment par une augmentation significative du moment 
magnétique, comme observé dans le cas de la monocouche Mol; dopée au Mn, où l'amélioration 
de ce paramètre est particulièrement marquée. Ainsi, le dopage substitutionnel apparaît comme 
une stratégie prometteuse pour améliorer les propriétés magnétiques des matériaux à base de 


Mo, avec des implications potentiellement importantes pour les applications en électronique et 


en spintronique. 


La même concentration a été utilisée pour le dopage de la monocouche VCI3 avec les 
éléments TM 3d (TM = Sc, Ti, Cr, Mn et Fe). On a constaté que les monocouches dopées au Sc 
et au Ti présentent un comportement semiconducteur ferromagnétique, avec des gaps d’énergies 
indirectes dans le canal de spin-up, tandis que la monocouche dopée au Cr présente des gaps 
d’énergies directes dans les deux canaux de spin, créant ainsi une monocouche semiconductrice 
ferromagnétique à gaps directes. Plus remarquablement, la monocouche VCI; dopée au Mn 
présente un comportement semiconducteur spin-gapless (SGS), tandis que celle dopée au Fe se 
révèle être un semiconducteur ferromagnétique bipolaire (BFMS). En outre, d'importants 
moments magnétiques de 4,25 us, 4,5 ue et 4,75 up par maille unitaire sont respectivement 
observés dans les cas dopés au Cr, au Mn et au Fe. De manière exceptionnelle, un couplage 
ferromagnétique est également observé entre le dopant Fe et les atomes de CI voisins, ce qui 


renforce la stabilité ferromagnétique de la monocouche VCH. 


Ces conclusions tracent un tout nouveau chemin vers l'optimisation des caractéristiques 
spintroniques des matériaux 2D dotés de leur propre magnétisme intégré. De plus, notre 
recherche offre la possibilité d'explorer l'incorporation d'autres éléments des groupes TM 4d et 
Sd, ouvrant ainsi la porte à la création de dispositifs électroniques et spintroniques de pointe, 


affichant des performances exceptionnelles. 
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Lorsque l'on effectue des calculs électroniques avec le code VASP, il génère de nombreux 
fichiers de sortie. L'outil VASPKIT peut aider à post-traiter ces fichiers pour les convertir en 
fichiers de données. Ensuite, on peut utiliser le logiciel gnuplot pour créer des graphiques de 
visualisation à partir de ces fichiers de données, représentant des propriétés telles que la densité 
des états et les diagrammes de structure de bande. Ci-dessous, on trouve deux scripts gnuplot qui 
permettent de dessiner facilement les diagrammes de la densité des états et de la structure de 
bande de la monocouche Mol, en entrant simplement la commande suivante dans le terminal : 


"enuplot nom_du_ script". 


B.1. Script pour tracer la densité d’états 


set terminal pdf enh size 8,8 font 'Times-New-Roman,20' 
set output ‘"pristine dos.pdf" 

set encoding iso 8859 1 

set size 1,1 

set multiplot layout 2,1 columnsfirst 

# --GRAPH a (Total density of states of VC13 monolayer) 
set xr [-2:5:3:.5] 

set yr [-58:58] 

set key at 2.5,58 

set key font 'Times-New-Roman, 22' 

set key spacing 1.0 

set label at -2.88,59 "(a)" font 'Times-New-Roman, 28! 
set border lw 2 

set style fill solid 1.0 noborder 

set arrow from 0,-58 to 0,58 nohead dashtype '-' lw 1 
set arrow from -2.6,0 to 3,0 nohead 1lw 1 

set origin 0.0,0.5 

set size 1,0.5 

set tmargin 2.5 

set lmargin 8 

set ytics font 'Times-New-Roman, 22" 

unset xtics 

plot './tdos.dat' u 1:2 w 1 lw 2 1c rgb "grey40" title 'Total','./tdos.dat' u 
1:3 w 1 lw 2 lc rgb ‘"grey40" notitle, 

unset label 

unset key 

unset ylabel 

unset tmargin 

unset lmargin 

unset arrow 

# GRAPH a (V-3d) 

set xr [-2.5:3.5] 

set yr [-58:58] 

set key at 2.5,48 

set key font 'Times-New-Roman, 22' 

set key spacing 1.0 

set origin 0.0,0.5 

set size 1,0.5 
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set arrow from 0,-58 to 0,58 nohead dashtype '-' lw 1 

set arrow from -2.6,0 to 3,0 nohead lw 1 

set border lw 2 

set ylabel 'DOS (states/eV)' font 'Times-New-Roman, 28' offset -4,-8 

set tmargin 2.5 

set bmargin 0.5 

set lmargin 8 

unset xtics 

unset ytics 

unset xlabel 

PIOEs "./PDOS V UP.dat' u 1:($6+$7+$8+$9+$10) w 1 lw 2 1e rgb "red" title 'V- 

saur ./BDOS Va DW.dat' u 1:($6+$7+$8+$9+$10) w 1 lw 2 1c rgb "red" notitle 

unset arrow 

unset label 

unset tmargin 

unset lmargin 

ii comes GRAPH a (C1-3p) 

set xr [-=2°5:3;5] 

set yr [-58:58] 

set key at 2.5,38 

set key font 'Times-New-Roman, 22' 

set key spacing 1.0 

set origin 0.0,0.5 

set size 1,0.5 

set arrow from 0,-58 to 0,58 nohead dashtype '-' lw 1 

set arrow from -2.6,0 to 3,0 nohead lw 1 

set border lw 2 

set tmargin 2.5 

set lmargin 8 

unset ytics 

unset ylabel 

plot './PDOS C1 UP.dat' u 1:($3+$4+$5) w L lw 2 le rgb ‘green" title 'Cl-3p', 

"./PDOS C1 DW.dat' u 1:($3+$4+$5) w 1 lw 2 le rgb ‘green" notitle 

unset arrow 

unset tmargin 

unset lmargin 

unset label 

# —--------- GRAPH b (V 3d states separated) 

set xlabel 'Energy(eV)' font 'Times-New-Roman, 28'offset 0,0 

set xt [=2-5:3;:5] 

set yr [-12:20] 

set key at 2.5,20 

set key maxcolumns 2 maxrows 3 

set key font 'Times-New-Roman, 22' 

set key spacing 1.4 

set origin 0.0,0.03 

set size 1,0.5 

set arrow from 0,-12 to 0,20 nohead dashtype '-' lw 1 

set arrow from -2.6,0 to 3,0 nohead lw 1 

set border lw 2 

set tmargin 0.6 

set lmargin 8 

set ytics font 'Times-New-Roman, 22" 

set xtics font 'Times-New-Roman, 22' 

unset ylabel 

PIOE "./PDOS V UP.dat' u 1:7 w 1 lw 2 1e rgb ‘magenta" title 'd x z/d y z' 
./PDOS V DW. dat' u 1:($7+$9) w L lw 2 le rgb ‘magenta" 
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necitlert -/PDOS-NOUP.datl u 1:6 w l'lw 2 1c rgb "blue" title "d x y/d {x2- 

y°}', './PDOS V DW.dat' u 1:($6+$10) w 1 lw 2 lc rgb "blue" 

noires ./BDOS VAUP dat ue 2 Cr brown title NiEz2i7 
./PDOS V DW.dat' u 1:8 w 1 lw 2 1c rgb "brown" notitle, 

unset arrow 

unset tmargin 

unset lmargin 

unset multiplot 


B.2. Script pour tracer la structure de bande 


set termopt enhanced 
set autoscale 
unset log 
unset label 
set xtic auto 
set ytic auto 
set key top right nobox 
set palette gray 
set encoding iso 8859 1 # Para poner acentos 
set grid 
set title ‘"" 
set xlabel "Pal" 
set ylabel "Energy (eV)" # font 
set key at 1.56,2.55 
set key font 'Times-New-Roman, 25' 
set key spacing 0.8 
set style line 1 LCL Iwe3 5 pt Ll'ps MIO  r9b red 
set style line 2 1t 1 lw 3.5 pt 1 ps 1.0 le rgb "web-blue" 
set style line 3 1t 1 lw 2.0 pt 1 ps 1.0 le rgb "black" 
set label at -0.1,3.5 "(b)'" font 'Times-New-Roman, 28! 
set arrow 1 nohead from 0.573,-2.5 to 0.573,3.5 dashtype 2 lw 1 
set arrow 2 nohead from 0.905,-2.5 to 0.905,3.5 dashtype 2 lw 1 
set arrow 3 nohead from 0,0.00 to 1.567, 0.00 dashtype 2 Llw 1.5 
set arrow from 0.95,0.29 to 1.1,0.75 lw 6 
# set arrow 4 from 0.0,-0.001137 to 0.0,2.355072 ls 3 
# set label at 0.05,1.9 "E {g}=2.356 eV" font 'Times-New-Roman, 28' 
set object circle at 0.905,0 size scr 0.033 fc rgb "black" 
set xr [0.0:1.567] 
set yr [-2.5:3.5] 
set xtics font "Symbol" 
set xtics ("{/Symbol G}" 0, "{M}" 0.573, "{/Symbol K}" 0.905, "{/Symbol G}" 
1.567) 
set style data histogram 
set style histogram cluster gap 1 
set style fill solid 
set term pdf size 8,8 enhanced font ‘Times New Roman, 28" 
set output "pure bs.pdf" 
set multiplot 
plot for [i=2:2] 'BAND.dat' using 1:(column (i)) title "up" with 1 ls 1, \ 
for [i=3:3] 'BAND.dat' using 1:(column (i)) title "down" with 1 1s 2, \ 


Permite pone “super/ {sub} indices 
scale axes automatically 

remove any log-scaling 

remove any previous labels 

set xtics automatically 

set ytics automatically 


RH + + 


unset xlabel 


Code source Gnuplot 


unset ylabel 

unset xtics 

unset arrow 

unset object 

set object 1 rect from graph 0, graph 0 to graph 1, graph 1 back 

set object 1 rect fc rgb 'white' fillstyle solid 1.0 

xsize= 0.5 # Controls the image x size in the canvas 

ysize= 0.25 # Controls the image y size in the canvas 
set size xsize, ysize 
set arrow 2 nohead from 0.905,-0.1 to 0.905,0.1 dashtype 2 lw 1 
set arrow 3 nohead from 0.8,0.00 to 1.0, 0.00 dashtype 2 Llw 1.5 
set xr [0.8:1.0] 
set yr [-0.1:0.1] 

set origin 0.46,0.55 

set object circle at 0.905,0.008 size scr 0.05 


plot for [i=2:2] 'BAND.dat' using 1l:(column (i)) title "" with 1 ls 1; 


unset style line 
reset 
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Aperçu sur le code de calcul 


C.1. Caractéristiques générales du VASP 


Tous les calculs rapportés dans cette thèse ont été effectués en utilisant Vienna Ab initio 


Simulation Package, plus communément connu sous le nom de VASP. VASP, initialement 


développé par Kresse et Hafner en 1994 et ensuite amélioré par Kresse et Furthmüller en 1996, 


est un logiciel écrit principalement en Fortran conçu pour effectuer des calculs quantiques ab 


initio [1,2]. VASP utilise soit des pseudopotentiels de Vanderbilt [3], soit la méthode des ondes 


augmentées par projecteur [4], en combinaison avec une base d'ondes planes. Voici quelques- 


unes des caractéristiques clés de VASP : 


V4 


> 


Effectue les minimisations électroniques et ioniques (relaxations ou optimisations) et 
calcule les propriétés de l'état fondamental (énergies totales, structures de bande, densité 
des états, spectres phononiques, etc.) en utilisant les orbitales KS étendues avec 
l'ensemble de bases d’ondes planes. 

Résout les équations KS de manière auto-cohérente avec une diagonalisation matricielle 
itérative et un mélange des densités de charge. 

Calcule les termes d'énergie de l'hamiltonien KS et les densités de charge dans les espaces 
réel et réciproque via des transformations de Fourier rapides (FFT). 

Fonctionne pour n'importe quel système, y compris les atomes, les molécules, les solides 
en volume, les surfaces, les amas, etc. 

Fournit la plupart des pseudopotentiels, y compris les pseudopotentiels ultra doux [3], les 
potentiels PAW [4], et ainsi de suite, pour tous les éléments. 

Fournit les fonctionnelles d’échange et de corrélation de LDA, PW91 [5], PBE [6], 
hybrides, etc. 

Fonctionne en mode parallèle (efficacement jusqu'à 32 nœuds) et peut traiter jusqu'à 4000 
électrons de valence. 

Fournit des dizaines de paramètres (tags) avec des valeurs par défaut et convient donc 


parfaitement à la fois pour un débutant en DFT et pour un expert. 


Pour exécuter VASP, plusieurs fichiers d'entrée essentiels sont nécessaires : 


INCAR : Ce fichier contient divers paramètres de contrôle pour le calcul, tels que les 


méthodes de relaxation électronique, les critères de convergence et d'autres paramètres. 
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POSCAR : Ce fichier définit la structure cristalline du système à étudier, y compris les 
positions des atomes, les paramètres de réseau et le nombre d'atomes dans la maille 
unitaire. 

POTCAR : Le fichier POTCAR contient essentiellement le pseudopotentiel pour chaque 
espèce atomique utilisée dans le calcul. 

KPOINTS : Ce fichier détermine la grille de points k pour l'échantillonnage de la zone de 
Brillouin dans l'espace réciproque. Il influe sur la précision et l'efficacité des calculs. 
Fichiers d'entrée facultatifs : Selon le calcul spécifique, vous pourriez avoir des fichiers 
d'entrée supplémentaires, comme WAVECAR (pour les reprises), CHGCAR (pour 


spécifier la densité de charge initiale) ou d'autres fichiers spécifiques à la tâche. 


C.2. Déroulement de VASP 


La Figure C.1 montre le déroulement général d'une exécution de VASP, comprenant à la 


fois les minimisations électroniques et ioniques avec les fichiers d'entrée et de sortie 


correspondants : 


V4 


Sélectionner un pseudopotentiel et une fonctionnelle d'échange-corrélation, puis 
construire une densité électronique d'essai. 

Utiliser une base d'onde plane pour l'expansion des orbitales KS, calculer l'énergie 
cinétique et le potentiel électrostatique externe V,,-(r) aux positions atomiques données 
et à la structure cristalline. 

Calculer l'énergie d'Hartree VA(r) et l'énergie d'échange-corrélation V,.(r), puis 
construire l'hamiltonien KS, Âxs. 

Résoudre les équations KS (une équation aux valeurs propres matricielle) par une 
diagonalisation itérative dans les espaces réel et réciproque. 

Avec les orbitales nouvellement calculées, générer une nouvelle densité électronique et 
répéter le processus jusqu'à l'atteinte de l'autocohérence. 

Une fois que la minimisation électronique est atteinte, calculer les forces sur les atomes et 
déplacer les atomes vers de nouvelles positions présentant des forces plus faibles. 

Répéter les minimisations électroniques et ioniques en série jusqu'à ce qu'un critère de 


convergence énergétique soit atteint. 
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INCAR, KPOINTS, POSCAR, POTCAR 


> Evaluer Îxs 


Résoudre les équations KS 
Diagonalisations itératives/sous-espaces 
Solution: Nouveau {W;(r)} 


| 


Nouvelle n(?) et F[n(r)] 


Convergé/No 
Convergé/Oui 


Relaxation ionique 


Convergé/Non 
Convergé/Oui 


CHGCAR, CONTCAR, WAVECAR, et ainsi de suite. 


Figure C.1 : Organigramme de minimisation 
électronique/ionique dans le code VASP. 
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